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ABSTRACT 

A  general  method  is  presented  for  finding  asymptotic  solutions  of 
initial-boundary  value  problems  for  linear  hyperbolic  partial  differential 
equations.   A  large  parameter  X   appears  in  the  equation,  multiplying  a 
lower  order  (dispersive)  term.   It  also  may  appear  in  the  initial  data 
and  the  inhomogeneous  (source)  term,  in  a  variety  of  ways.   This  gives 
rise  to  a  variety  of  different  types  of  asymptotic  solutions.   The  expan- 
sion procedure  is  a  "ray  method",  i.e.,  all  of  the  functions  that  appear 
in  the  expansion  satisfy  ordinary  differential  equations  along  certain 
space-time  curves  called  rays.   These  rays  do  not  lie  on  characteristic 
hypersurfaces,  but  instead  fill  out  the  interior  of  the  characteristic 
hypercone.   They  are  associated  with  an  appropriately  defined  group  veloc- 
ity.  The  details  of  the  expansion  are  presented  for  a  simple  second  order 
hyperbolic  equation.   The  applicability  of  the  method .to  symmetric  hyper- 
bolic systems  and  to  the  integro-differential  equations  for  dispersive 
dielectrics  is  discussed. 


1.   Introduction 

In  recent  years  asymptotic  methods  have  been  developed  for  the  solution 
of  certain  boundary-^^alue  problems  and  initial-value  problems  for  partial 
differential  equations.   These  problems  involve  a  parameter^  and  the  methods 
provide  one  or  more  terms  of  the  asymptotic  expansion^  say  for  large  values 
of  the  parameter  J  of  the  solution  of  the  problem.   They  are  often  applicable 
to  problems  for  which  exact  solutions  are  not  known ^  and  even  for  problems 
which  can  be  solved  exactly  it  frequently  happens  that  only  the  asymptotic 
expansion  of  the  solution  is  sufficiently  simple  to  be  useful  in  practical 
applications.   Furthermore ;,  it  is  invariably  true  that  the  methods  which 
yield  the  asymptotic  expansion  directly  are  very  much  simpler  than  the  pro- 
cedure which  involves  first  finding  the  exact  solution,  and  then  its  asymp- 
totic expansion. 

An  important  class  of  asymptotic  methods  is  characterized  by  the  fact 
that  certain  curves _,  often  called  "rays",  play  a  central  role  in  the  theory. 
The  rays  are  of  fundamental  importance  because  all  of  the  fixnctions  which 
make  up  the  various  terms  of  the  expansion  can  be  shown  to  satisfy  ordinary 
differential  equations  along  these  curves.   Thus,  in  a  sense,  the  method  is 
one  which  reduces  partial  differential  equations  to  ordinary  differential 
equations. 

The  "ray  method"  has  been  extensively  developed,  primarily  by 
J.  B.  Keller  and  his  co-workers  at  New  York  University,  for  the  reduced 
wave  equation 

^  V  +  k^n^(X)v  =  0  (1) 

and  related  elliptic  equations.  When  applied  to  (l)  the  method  yields  a 
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"geometrical  theory  of  diffraction"'-  -'  which  generalizes  the  classical  theory 
of  geometrical  optics.   The  extensive  literatiire  on  the  asymptotic  theory  of 
the  reduced  wave  equation  and  Maxwell's  equations  has  recently  been  unified 
and  summarized  in  [2].   The  reduced  wave  equation  is,  of  course,  obtained 
from  the  wave  equation 

Au-     1 


c2(X)     **  (2) 

by  the  substitution  u(t,X)  =  v(X)e""^   .   Then  k  =<^/c^  and  n(X)  =  c^/c(X). 

(c  is  a  convenient  constant.   The  large  parameter  is  u  or  k.)   Similarly 
o 

the  "symmetric  hyperbolic"  system  of  equations 

n 
H  A^-|-.Bu=0,    (x^  =  t),  (5) 

v=0        V 
is  reduced  by  the  substitution  u(t,X)  =  v(X)e"      to  the  system 


H  A^-|^  +  (B-iu)  )v=0.  {k) 

V=l        V 

Here  X  =  (x  ,  . . .  ,  x  ) ,  A  and  B  are  real  kXk  matrices  which  are  functions 
of  X,  A°  is  positive  definite,  and  A  ,  . . . ,  a'^  are  symmetric.   The  asymptotic 
theory  of  {k)    for  large  w  has  been  studied  in  [J]. 

The  original  hyperbolic  system  (5)  contains  no  parameter.   However, 
P.  D.  Lax  '-  ^    has  obtained  the  asymptotic  solution  of  the  initial-value 
problem  for  (j)  with  "oscillatory  initial  data"  of  the  form 


u(0,X)  =  e^^^^^)  z(X).  (5) 

In  Lax's  problem  the  large  parameter  \  appears  only  in  the  initial  data, 


and  specifically  in  the  form  (5).   The  role  of  the  rays  is  played  "by  the 
"bicharacteristics"  which  are  space-time  curves  that  generate  the  character- 
istic hypersurfaces  of  (j). 

In  this  paper  we  shall  examine  a  method  which  is  applicable  to  the 
solution  of  initial  value  and  initial-boundary  value  problems  for  sym- 
metric hyperbolic  systems  of  the  form 

n 

zIa^I^  +  XBu=  f(t,X;x)  (6) 

v=0      V 

with  initial  conditions 

u(0,  X)  =  u^(Xj  X).  (7) 

The  large  parameter  X  appears  in  the  system  of  differential  equations 
(6) ,   multiplying  the  undifferentiated  term  Bu.   It  is  well  known  (See  [5] , 
p.  191  )  that  the  presence  of  such  terms  leads  to  the  phenomenon  of  dispersion*^ 
and  the  large  paremeter  in  the  dispersive  term  ensures  that  these  phenomena 
appear  in  the  asymptotic  solution.   For  this  reason  the  asymptotic  method 
when  applied  to  (6)  leads  to  results  which  are  quite  different  from  the 
solution  of  Lax's  problem.   For  the  dispersive  system  (6)  the  bicharacter- 
istics no  longer  serve  as  "rays."  In  fact  the  rays  do  not  lie  on  character- 
istic hypersurfaces,  and  they  represent  trajectories  of  points  which  move 
with  an  appropriately  defined  "group  velocity."  Thus,  for  example,  the 
asymptotic  construction  of  the  Riemann  function,  which  is  defined  by  (6)  and 

(7)  with  f  =  0  and  u  =  6(X) ,  leads  to  ravs  that  emanate  from  the  origin  in 
o 

space-time  and  fill  out  the  interior  of  the  characteristic  conoid.   The 
parameter  X  also  appears  in  the  source  term  f  of  (6)  and  the  initial  data 


For  dispersive  systems,  different  frequency  components  of  a  wave  travel  with 
different  velocities;  this  leads  to  distortion  of  wave  forms  and  other  inter- 
esting physical  phenomena. 
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function  u  of  (7).   ¥e  will  see  that  the  simplest  type  of  expansion  is  ob- 
tained when  u  is  oscillatory,  i.e.,  of  the  fomi  (5).   However  it  is  possible 
to  obtain  asymptotic  solutions  of  (6,7)   when  X  appears  in  u  and  f  in  a  variety 
of  different  ways.   This  leads  to  a  variety  of  different  types  of  asymptotic 
solution,  each  of  which  is  characterized  geometrically  by  a  typical  ray 
pattern.   Some  of  these  will  be  examined  in  sections  3  and  h   of  this  paper. 

An  interesting  example  of  a  dispersive  symmetric  hyperbolic  system  of  the 
form  (6)  is  given  by 

VXH-6^=J,  VXE  +  |i|^=0,  ^  -   e  f  ^  =  0.  (8). 

The  first  two  equations  are  simply  Maxwell's  equations  for  the  electric  and 
magnetic  fields  E,  H  in  M.K.S.  units.   The  third  equation  is  obtained  by 
assuming  that  the  current  J  is  entirely  due  to  the  motion  of  electrons  with 
charge  e,  mass  m  and  number  density  t)(X)  ,  and  that  these  electrons  satisfy 
Newton's  equation  with  the  force  term  given  by  eE.   In  (8)  (j)  is  the  "plasma 
frequency"  defined  by 

4'  =  ^  •  (9) 


The  system  (8)  describes  the  propagation  of  electromagnetic  waves  in  certain 
plasmas  such  as  the  ionosphere.   For  the  case  of  harmonic  time  dependence 
(8)  reduces  to  the  well  known  equations  of  the  (isotropic)  magneto-ionic 

theory'-  -' .   Let  us  introduce  the  velocity  of  light  c  -(fe|J.)  ^,   an  average 

I  * 

plasma  frequency  (p  ,  and  a  characteristic  length  a.   We  then  replace 

E,H,J,X,t  in  (8)  by  E',  H',  J',  X',  t'  and  set 
]/e'  E'  =  E,Vi^  H'  =  H,  J'='Ve~  (j)  J,  X'=aX,  ct'=at,  <j)(X)=(f)  b(X).       (lO) 


For  problems  with  boundaries  a  denotes  a  typical  boundary  dimension.   For 
problems  with  non-constant  ^{xj   a  is  an  average  value  of  the  function 

1^1 /'^'v  I    ,   Thus  (f)(X)  changes  by  a  small  fraction  of  itself  over  distances 
'  '    '      small  compared  to  a. 
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Then  (8)  becomes 


E   -  V  X  H  +  Xb(X)j  =0,   H^  +  V  X  E  =  0,  J^  -  Xb(X)E  =  0. 


(11) 


The  (large)  parameter  X  is  given  by 


X  =  -^  .  (12) 

X  and  t  are  dimensionless  space  and  time  variables  and  b(X)  is  a  dimensionless 
function  proportional  to  the  plasma  frequency. 

It  is  now  an  easy  matter  to  verify  that  (ll)  is  of  the  form  (6)  with 
f  =  0.   In  this  case  u  is  the  column  vector  with  components 


{■S,^,Y.^,Y.y^^,n^,^yJ^,J^,J^)   , 


(13) 


A°  is  the  9  X  9  identity  matrix,  and  A"'", A  ,P?   are  9  X  9  matrices  defined  by 


1 


k  A^ 


'O  -(K)   0 

(K)    0    0 

0    0    0 


0  0  bl, 
0  0  0 
-bl  0    0 


H^) 


Here  I  is  the  5  X  J  identity  matrix  and  (K)  is  the  3  X  J  matrix 


(k)  = 


0 

-S 

^^2 

^ 

0 

-^1 

[-^■2 

\ 

0  _ 

(15) 


Since    (k)   is  anti-symetric  we  note  that  A     is   symmetric  and  B  is  antisymmetric. 
The  last  condition  implies  that  the   system  conserves  energy  and  is  important 
for  our  expansion  procedure  as  we  will  soon  see. 

H  and  J  are  easily  eliminated  from   (ll)  to  obtain  the   second  order 
system  of  equations 
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E,,-  -^ 


+  V(V.E)  +xV(X)E  =  0.  (16: 


If  b  is  independent  of  x  and  E  =  li{x-.,x^)A,   where  M-  is  a  scalar  function 
.A  is  a  unit  vec 
nd  (l6)  reduces  to 


and  A  is  a  unit  vector  in  the  direction  of  the  3  axis,  then  V.  e  =  0 


u^^-  c^(X)Ai  +xV(X)u  =  0,  (17) 

with  c  =  1. 

Equation  (17)  is  a  simple  example  of  a  dispersive  hyperbolic  equation 
that  conserves  energy,  and  it  serves  admrably  as  a  model  for  demonstrating 
our  method  of  asymptotic  solution.   In  Section  (2), (3),  and  (k)   we  shall 
develop  methods  of  obtaining  asymptotic  solutions  of  initial  and  initial- 
boundary  value  problems  for  (I7)  and  for  the  same  equation  with  an  inhomogeneous 
(source)  term.   In  our  work  we  shall  deal  only  with  (17)  and  in  sections 
(3)  and  (k)   we  shall,  for  simplicity, further  restrict  oi;ir  considerations 
to  the  case  of  one  space  dimension.   However,  we  emphasize  that  the  methods 
are  quite  general.   Research  currently  in  progress  has  established  that  they 
are  applicable  to  symmetric  hyperbolic  systems  of  the  form  (6)  (with  anti- 
symmetric B)  in  any  number  of  dimensions.   The  main  condition  necessary  for 
application  of  the  method  can  best  be  described  in  terms  of  the  "dispersion 
equation" 

5(k^,  ...,  k^,  w)  =0  (18) 

which  can  be  obtained  by  requiring  that  a  plane  wave  of  the  form 


n  ^ 

u  =  const,    exp  ^    i>^(  V"     k„x  -wt)  )  (I9) 


expl    i^(  5~     ^v^v-^)  ( 


v=i 

"be  a   solution  of  the  partial  differential  equation  with  constant  coefficients. 
Then  the  condition  is  that  the   solutions 

tJ  =  h(k^,...,k^)  (20) 

of  (18)  be  real  for  real  k  , . . . ,k  .    (To  what  extent  this  requirement  can 
be  relaxed  is  not  yet  clear.  ) 

For  (17) 


-VcV-ft^    .^2^^  ^2  .  ^21; 


and  for  (6)  (with  constant  coefficients  and  f  =  O)  the  dispersion  equation 
takes  the  form 

n      ^ 
det  G  =  0;    G  =  Q-A  w  ;  Q  =  y  k/  -iB.  (22) 

v=l 

That  this  equation  has  only  real  solutions  co  for  real  k  follows  from  the 
fact  that  A  is  positive  definite  and  Q  is  hermitean^  for  w  is  an  eigenvalue 
of  the  hermitean  matrix  (A°)"  '^   Q(A  )"  '  . 

Before  proceeding  to  the  detailed  discussion  of  our  method  we 
would  like  to  point  out  that  it  is  also  applicable  to  the  very  important 
case  of  dispersion  of  electromagnetic  waves  in  dielectrics.   This  case  is 
not  included  in  (6)  because  the  electric  and  magnetic  field  vectors  E,  H 
in  a  dispersive  dielectric  satisfy  the  system  of  differential  and  integral 
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* 

equations 


VXE+i|S     =0,    VxH-i     ^  =0, 
c    Ot  c      at  ' 


i 


D(t,x)    =  E(t,X)    +  /      f(T,X)E(t-T,X)dT.  (23) 


Here  we  have,  for  simplicity,  taken  the  magnetic  permeability  to  be  U-  =  1. 
The  dielectric  permeability  e  is  a  function  of  frequency,  and  f  is  given  by 

f(t.X)=~y^    e-^"^[e{^,X)-l]d.^.  (2k) 

Since  e  is  dimensionless  it  must  be  a  function  of  ^/'^^^       where  w  is  a  characteristic 
frequency  of  the  dielectric.   We  now  replace 

X,  t,  w^e^f      by     X*,  t',  W,  e«,  f 
and  introduce  a  characteristic  length  a.   Then  we  set 

X'=  ax,    ct'    =  at,   0)'    =  w^^,    e'(^',V)   =e(w,X),   af '  (t ' ,  T  )  =  cf  f  t  ,X) . 

(25) 

With  this  transformation   (23)  and   (2^1-)  become 

00 

VxE+^  =  0        ,      VxH-^=0,      D(t,x)    =  E(t,x)    +  /  f(T,X)E(t-T,X)dT, 


P  (26) 

f(t.X)   =^    J       e-^"fe(a.,X)_i]dc..  (27) 


X   =  -2-      .  (28) 


c 


*  See    [b],    section  58. 


Although  the  system  of  equations  (26)  is  quite  different  from  (6)  a  very 
similar  asymptotic  method  is  applicable  to  the  solution  of  (26)  for 

■X- 

large  ^. 

Let  us  indicate  briefly  the  steps  involved  in  the  asymptotic  method. 
For  problems  which  can  be  solved  exactly   ,  examination  of  the  asymptotic 
expansion  of  the  solution  shows  that  it  consists  of  a  sum  of  terms,  each  of 
which  is  an  asymptotic  power  series  in  X'   involving  a  "phase  function"  and 
an  infinite  sequence  of  "amplitude  functions".   For  complex  problems,  we 
therefore  assume  that  the  solution  is  also  given  by  a  sum  of  such  series. 
By  inserting  such  a  series  into  the  partial  differential  equation,  we  find 

first  that  the  phase  function  satisfies  a  first  order  partial  differential 

It*** 
equation  which  can  be  solved  by  the  method  of  characteristics   ,   The 

characteristic  curves  are  of  course  the  rays  which  we  have  mentioned.   They 

are  space-time  curves  .   The  characteristic  (ordinary)  differential  equations 

will  be  called  "ray  equations".   Thus,  by  the  method  of  characteristics  the 

phase  function  can  be  obtained  by  integrating  an  ordinary  differential 

equation  along  the  rays.   We  also  find  that  the  amplitude  functions  satisfy 

ordinary  differential  equations  along  the  rays  and  that  these  "transport 

equations"  can  be  solved.   In  order  to  find  the  rays,  and  the  phase  and 


*  We  hope  to  complete  our  research  on  systems  of  the  form  (6)  and  (26)  in 
the  near  future. 

**-  For  example,  problems  with  constant  coefficients  and  no  boundaries. 
***See  [5],  Chapter  II. 


amplitude  functions,  it  is  necessary  to  specify  initial  conditions  for  all 
of  these  ordinary  differential  equations.   In  some  cases  (for  example,  the 
case  of  oscillatory  initial  data)  the  required  initial  conditions  follow 
directly  from  the  data  of  the  prolDlem.   In  others  the  initial  conditions 
are  obtained  from  the  solution  of  a  "canonical  problem".   A  canonical 
problem  is  one  with  the  same  local  features  as  the  given  problem.   It  is, 
however,  sufficiently  simple  to  be  solved  exactly.   The  required  initial 
conditions  for  the  given  problem  are  obtained  by  examination  of  the 
asymptotic  expansion  of  the  solution  of  the  canonical  problem. 

It  is  by  now  clear  that  the  asymptotic  method  involves  several 
unproved  assumptions.   It  is  therefore  reasonable  to  ask  whether  it  can  be 
proved  that  it  does  Indeed  yield  the  asymptotic  expansion  of  the  exact 
solution  of  the  given  problem.   No  complete  proof  of  this  fact  has  been 
given.   Nevertheless  there  is  abundant  evidence  of  the  validity  of  the  method. 
This  evidence  has,  so  far,  been  obtained  in  two  ways: 

1.   Rroblems  for  (17)  in  one  space  dimension  can  be  solved  exactly 
by  fourier  analysis  with  respect  to  t.   The  resulting  ordinary  differential 
equation  (with  independent  variable  x)  can  be  solved  asymptotically  by 
standard  methods  for  ordinary  differential  equations  (W.K.B.  method).   Then 
the  asymptotic  expansion  of  the  resulting  Integral  representation  can  be 
obtained  by  the  "method  of  stationary  phase"  and  related  methods.   In  every 
case  that  has  been  examined  the  resulting  asymptotic  expansion  of  the  exact 
solution  agrees  perfectly  with  the  results  obtained  in  this  paper  by  the 
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"ray  method".   These  investigations  yield  not  only  confidence  in  the  validity 
of  the  ray  method  but  also  an  appreciation  of  the  great  simplicity  of  that 
method  when  compared  with  the  approach  just  outlined. 

2.   In  the  case  of  initial  value  problems  for  (6)  with  oscillatory 
initial  data  it  is  possible  to  estimate   the  remainder  of  the  asymptotic 
series  by  energy  inequalities.   This  approach  is  similar  to  one  outlined  by 
Lax  for  the  initial  value  problem  (3j5);  and  yields  an  indirect  proof  of 
the  asymptotic  nature  of  the  expansion  in  this  case. 

The  ideas  developed  in  this  paper  had  their  origins  mainly  in  the 
work  of  J. B.  Keller  [l,2]  and  G. B.  Whitham  [6] .   The  latter  approached  the 
problem  of  the  propagation  of  dispersive  waves  from  a  physical  point  of  view, 
with  considerable  success.   Our  results  in  sections  (3.^)  and  (^.1)  are  easily 
seen  to  be  equivalent  to  Whitham' s  results.   In  attempting  to  develop  Whitham' s 
approach  into  a  systematic  expansion  procediore,  we  first  introduced  the  large 
parameter  ^   as  in  (6),    (ll),  (l?)  and  (26,27).   This  enabled  us  to  proceed  in 
a  manner  very  similar  to  Keller's  method  for  the  reduced  wave  equation.   The 
essential  idea  of  the  "cannonical  problem"  also  had  its  origins  in  Keller's  work. 


*  The  author  is  indebted  to  Mr.  N.  Bleistein  who  has  carried  out  these 
comparative  studies  and  is  currently  engaged  in  other  critical  examinations 
of  the  ray  method. 

*^  So  far  we  have  been  able  to  obtain  estimates  in  the  Lp  norm(energy  norm). 
It  is  not  clear  whether  pointwise  estimates  can  be  obtained  by  these  methods. 
***  It  is  worth  noting  that  such  indirect  proofs  are  much  more  difficult,  if 
not  impossible  in  the  case  of  elliptic  equations  such  as  (l)  and  (k) . 
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Some  features  of  the  expansion  we  shall  employ  are  similar  to  those 
of  expansions  that  have  been  used  for  the  time- dependent  Schroedinger  equation. 
In  fact  the  mechanical  terms  which  we  shall  use  such  as  "hamiltonian", 
"lagrangian",  etc.  occur  quite  naturally  when  the  expansion  is  applied  to  the 
initial  value  problem  for  the  Schroedinger  equation.   In  this  case  the 
appropriate  large  dimensionless  variable  is  ^  =  r*-  where  a  is  a  characteristic 
length  p  is  a  characteristic  momentum  and  h  is  Planck's  constant.   We  have 
found  it  convenient  as  well  as  suggestive  to  continue  to  use  the  mechanical 
terms  for  all  our  equations. 

Before  proceeding  to  the  asymptotic  solution  of  (17)  we  add  one 
final  remark:   If  we  introduce  the  transformation 

x«  =  Xx,   t'  =  U,  u'  (f  ,  x'  )  =  u(t,x)  (29) 

(17)  becomes 

u'^^,-  [c'(x»)]^  A^,u'  +  [b'(x')]^u'  =0.  (30) 

Here 

c'(x')  =  c(x)  =  c(x'A),  b'(x')  =  b(x)  =  b(x'A).  (31) 

Thus,  in  finding  asymptotic  solutions  for  X'^<=°    of  problems  for  (YJ),   the 
transformation  (29)  shows  that  we  will  also  be  obtaining  asymptotic  solutions 
for  large  t',  x'  of  corresponding  problems  for  equation  (30)  in  which  no 
parameter  appears.   It  should  be  pointed  out,  however,  that  (3I)  implies  that 


.1^- 


^  =o(i),  ^  =oih.  (32) 

Therefore^  for  the  case  of  variable  coefficient  functions,  the  asymptotic 
solutions,  as  applied  to  (30)  are  valid  only  for  slowly- varying  coefficient 
functions.   If  boundaries  are  present  in  the  problem  for  (17)  then  (29) 
shows  that  the  corresponding  boundaries  for  (30)  must  have  large  dimensions. 
In  all  of  the  succeeding  sections,  vectors  are  denoted  by  capital 
letters.   We  also  employ  the  surrmiation  convention  with  respect  to  repeated 
indices. 

jy2        Asymptotic  Series  Solutions  of  a  Hyperbolic  Equation 

Let  us  assume  that  a  solution  u(t,X)  =  u(t,x^,. . .x^)  of  (1.15)  has 

an  asymptotic  expansion  of  the  form 


u  ~  e 


.As(t,X)  y-   (.^)-m  ,(tn)(^^j^j.  Jm)  ^  q  f or  m  =  -1,-2,...  .         (l) 


By  inserting  (l)  in  (I.I5)  and  equating  to  zero  coefficients  of  (lA)"   we 

■it- 
obtain  the  dispersion  equation 

o^{X)s       s       -  s  s   +  b2(x)  =  0,  (2) 

and  the  transport  equations 


2(s   z 


Here   Qs  =  s     -  c'^  s 

V  V 


(^)  .  c-%z^^)  ;  z(™^  Ds  +  Q  z(^-l^=  0;  m=0,l,2,  ....     (3) 


*Here  the  dispersion  equation  is  a  partial  differential  Equation  for  s  whereas 
in  section  1  it  was  an  algebraic  equation.   No  confusion  will  arise  from  the 
double  meaning  of  the  term. 
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We  choose   a  constant  reference  spaed  c  and  introduce  the  following 
notation,  which  will  be  used  throughout  the  paper: 

'^o  =  ^o*'  S  =  %  '  ^o  =  \  =  -  f; '  ^  =  -%'  ^  =  (ki,...,y, 

K  =  (k^,K),  X  =  (x^,X),  k^  =  k^k^,  ^^  =  ^o  ■"  ^^• 

Then  the  dispersion  equation  takes  the  form 

c\^-  </  +-b^  =   0,  (5) 

or,  if  we  solve  for  u, 

{J  =  h(X,K)  =  i-  h^(X,K);  h^(X,K)  =  Vc^^*^.         (6) 

For  each  of  the  two  values  of  the  hamiltonianj  h,  the  dispersion  equation  can 
be  written  in  the  form 

h(X,K)-to  =  h(X,K)-»c^k^  =  0.  (7) 

We  now  define  the  quantities 
8v  =  ^  =  "IT-  '  So  =  V  ^=   (S^f'-'>S^^)>     G  =  (g^,  G), 


/v2    2     2    2  2,^  ^2,  2. 

g  =  Sq  +  g  *   g  =  gygy  =  C  (1-h  /u  ). 


(8) 


G  is  the  group  velocity  vector,  and  g  is  the  group  speed.  By  analogs^  with  plane 
wave  solutions  of  (l.l?)  we     see  that  it  is  appropriate  to  refer  to  w  and  K 
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as  (angular)  frequency  and  propagation  vector,  or  if  we  wish  to  emphasize 
the  mechanical  analogy  we  may  refer  to  K  as  the  moment im  vector. 

In  the  following  it  will  be  noted  that  we  shall  make  no  essential  use 
of  the  explicit  form  (6)  of  the  hamiltonian  function.   For  this  reason  it 
is  not  surprising  that  o\:ir  methods  can  be  applied  to  a  large  class  of  partial 
differential  equations,  as  we  have  asserted  in  the  introduction. 
2.1  the  ray  equations  The  first  order  partial  differential  equation  (Hamlltc 
Jacobi  equation)  (t)  can  now  be  solved  by  introducing  the  characteristic 
equations  (Hamilton's  equations  of  motion) 

(9) 


*  =  ...  -  =  =,i,  =  0,  .,  -  -  1^ 


Since  t  =  1,  we  may  choose  the  independent  variable  in  this  system  of  ordinary 
differential  equations  to  be  t.   Each  solution  of  (9)  defines  a  curve  X(t) 
in  X- space  or  equivalently  a  space-time  curve  [t,x(t)]  which  we  shall  call 
a  ray.   Prom  (h)   and  (9)  we  now  have 


g  =  s   +  s   k^=  £  (10) 

dt    t    X,,  ^ 


where  \    is  the  lagrangian, 


v^--  (^^^ 
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2,2  solution  of  the  dispersion  equation 

An  Initial  value  problem  for  s(t,X)  can  now  be  solved  by  means  of 


(9)  and  (n).  Let  Z   =  {a^, . . . ,    o  ),   and  let 


t  =  t^(E),  x^  =  x^^(^.);  V  =  1,  ...,  n  (12) 


be  the  parametric  representation  of  a  hypersiirface  on  which  s  has  the  prescribed 
values 

s[t,(t),   x^„(E)]  =  s^(E).  (13) 

Then  differentiation  of  (13 )  with  respect  to  a .   yields 

J 

3x,^      N,o    Ss 

■s-337  -"%  =s5^;o=i,  ......  m 

(l4)  and  (6)  provide  ji+1  algebraic  equations  for  the  determination  of  u  =  to  (e) 
and  k  =  k  (E),  v  =  1,  ...,|a,  at  each  point  L        of  the  initial  hypersurface 
(12).  Then  t  (E),  x  (e),co  (I.),  k  (E),   provide  initial  values  for  the  solution 
of  (9).  If  this  solution  is  de.iiJted  by  x  =  x  (t;  S),  k  =  k  (t;  E), 
c  =  u)(t;E)  =  w(E)3  then  x^[t^(E);E]  =  x^(e),  \[t^(Z)f^]    =  ^vo^^^'  ^^"^ 
w(E)  =  u  (2),  In  terms  of  this  solution  (lO)  may  be  integrated  to  yield 


[t,x/t;E)]  =  s^(E)  +  f         >^[\(*S  ^)*  k^(tSz)Jdt«.        (15) 


Correeponding  to  the  two  values  of  (6),  (15)  yields  the  two  solutions  of  the 
initial  value  problem  for  s. 
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The  initial  value  problem  for  s  is  considerably  slnjilified  if  the  initial 
hypersurface  is  the  hyperplane  t  =  t  .  In  this  case  we  may  determine  k  (Z) 
and  w  (E)  at  once.  We  set 

t^(E)  =  t^,  x^(£)  =  cr^,  V  =  1,  ...,  n,  (16) 


in  (12).  Then  (ik)   and  (6)  yield 


It  is  interesting  to  note  that  (ll)  implies  that  s  is  constant  on 
rays  if  and  only  if   k  ^  =  h.   But  this  is  equivalent  to  the  condition  that 
h  he  homogeneous  of  degree  one  in  K.   In  this  case  the  original  partial 
differential  equation  contains  no  lower  order  terms  and  is  non- dispersive. 
For  non- dispersive  equations  it  is  easily  seen  that  the  dispersion  equation 
(as  a  differential  equation  for  s)  is  the  "characteristic  equation",  the 
level  surfaces  of  its  solution  s(t,X)  are  characteristic  hyper- surfaces  of 
the  original  equation,  and  the  rays  lie  on  the  characteristic  surfaces.   None 
of  these  statements  is  true  for  dispersive  equations  such  as  (l,17). 

2,3  solution  of  the  transport  equations 

V7e  turn  our  attention  now  to  the  transport  equations  (3)  which  may  be 

written  in  the  form  ,,0. 

(lo; 

If  we  set  z^   B  2,  the  zero  order* equation  is 


2  2 

z  +  (i7Tns)z  =  z  +  g  z     +  (|-  a  s)z  =  0.  (19) 
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These  are  ordinary  differential  equations  along  each  ray.   To  solve  the  homogenous 
equation  (I9)  we  first  note  that 


X^Xy       tt      t    dx^     *^0  c    ^     '^v  c 


where 


(20) 


v=(|j^,  V),  ^  .  (g^^  a),  g,=  c^  .  (21) 

Then  (19)  can  be  written  in  the  form 

2|  6^.  Oz  +z^.  (iii^^)  =0,  (22) 

c  c 

2 
2-G)  =0  (23) 

c 

For  an  arbitrary  domain  D  in  E-spa^e  we  now  ajrply  Gauss*  theorem  to 
the  "ray-tube"  In  30- space  determined  by:  E  In  D,  t,  <  t  <  tgt 


N 


Since  G  »  (c^j^CJ)  Is  everywhere  parallel  to  the  irays  X  «=  (x^X)  =  [c  t,  X(t,S)], 
(23)  yields 

2 


[t^2   t=*^   J 

(e  In  d)  (2:  m  D) 


0     ..  ..J 


Here  J  is  the  Jacoblan, 


a(x^,...,x  )                pec  (t,E)] 
j(t,E)  =    -^ B    =  det    -g^ .  (25) 


V""n 
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2. 

Since  D  was  arbitrary  and  w  is  constant  on  a  ray  (2k)   implies  that  ~^  is 

c 
constant  on  a  ray  (provided  j  does  not  vanish).   Thus  the  solution  of  (19)  is 

given  by 


tUfr  [m-} 


1/2 


^(t)  =  zct^)  ^>^  lWrr'^  (26) 


Here  z(t)  =  z[t,X(t,E)],  c(t)  -  c[X(t,E)],  and  j(t)  «  j(t,2).  The  solxrfclon  of 
the  Inhomogeneous  equation  (18)  can  now  "be  obtained  by  the  method  of  "variation  of 
parameters".  The  result  is 


(27) 

m  =  0,1,2,,,., 


2,k  Constant  coefficients 


The  foregoing  results  siniplify  considerably  in  case  the  coefficient 
functions  are  constant:  b  =  b  ,  c  =  c  ,  In  this  case  the  ray  equations  (9) 
imply  that  k  .u,g^,  and  ^  are  constant  on  each  rayj  i.e. 

(28) 

^=io(^)'^oSvo-V 

Purthermore  the  rays  are  the  straight  lines 


x^  »  x^(E)   +  g^(E)    p-t  JE)]    .  x^^-  t^g^   +  g^^t,  (29) 

aod 

8  =  8 Jr.)  +  iJj:)[t-tjE)]  =  8^-  t^  i^  -^     ^^t.  (30) 

(2  5)  and  (  29)  now  linply  that  (in  general) 

j(t,E)  =  det   ^(x^-  t^g.^)^^.  +  (g^^J^^tl  (31J 

is  a  polynoinial  of  degree  n  in  t,  vith  coefficients  which  are  functions  of  E. 
Hence  if  t,(E),,,,,t  (E)  are  the  roots  of  this  polynomial  we  may  write 

j(t,S)  =a(E)  •ff''  [t-T  (e)J,  (32) 

and  (27)  takes  the  special  form 


°*o 


J(t„)  (VTj.)....(W 

TrE7  -  w~r^)....(t-Tj}   •  (31.1 

2.^  plane  waves 

The  siii5)lest  exanrples  of  our  formal  expansion  are  the  pQane  waves. 
We  assume  that  the  coefficients  are  constant  and  talte  the  initial  surface  to 
he  the  hyperplane  t  =  0,     For  arbitrary  constant  K  »  (k^,...,k  )  we  set 
s®(E)  B  kyffy  .     fflien  (17)  yields 


We  define  g^  and  £^  by  ^8)  .  Then  (29),  yields  the  family  of  parallel 

(36) 


'vo 

straight  rays 


and  (30)  yields 


^  =  ^v%  +  (^vo^vo  -  "o^*  =  ^v^v  -  "o*-  (3TJ 


Prom  (36)  we  see  that  j  5  1.  We  now  set  z^^'   =  const,  on  the  plane  t  »  0  for 

m  =  0,1,2,..,.  Then  (33)  iugplles  that  z(^'(t,X)  Is  constant.  If  we  Insert  these 

results  in  (l)  we  obtain 

u(t,X)  =  2(X)e^'^(^V^v  -  ''o*^  zW  ^   2^(l7v)-^z("'\     (38) 

In  fact  (38)  is  an  exact  solution  of  (1.15)  with  arbitrary  z(?v)  provided  w 
satisfies  (35). 

2,6  waves 

Any  asyn5>totic  series  solution  of  (l.lT)  of  the  form  (l)  will  be 
called  a  wave.  Plane  waves  are  a  special  case.  We  have  seen  tteit  if  the  value 
of  s  is  prescribed  on  a  siiitable  ii-dimensional  manifold  M,  two  solutlcms  s  of  the 
dispersion  equation  are  determined  (at  least  in  a  neighborhood  of  M).  If  the 
values  Of  the  functions  z^°^'  for  m  «=  0,1,2,.,,  are  also  prescribed  on  M  exactly 
two  waves  are  coo^letely  determined.  In  sticceeding  sections  we  will  atteuqprt  to 
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show  how  a  svun  of  one  or  more  waves  can  be  constructed  to  yield  asymptotic 
solutions  of  suitable  problexne  for  (l,l7), 
2,7  lover  dimensional  initial  manifolds 

For  some  problems  it  viU  be  necessary  to  ccaistruct  waves  which  are 
determined  by  the  values  of  s  and  z^"^'  on  a  manifold  M  of  dimension  6  <  u. 
In  the  extreme  case,  o  =  0,  M  Is  a  single  point  (^q^X  )  «=  (t  ,  X-q,,,.,X  ), 
and  s  has  the  prescribed  value 

«CV^o)  =  ^o-  (39) 

We  shall  then  require  waves  for  which  all  of  the  assosciated  rays  pass 
through  M.   These  rays  are  obtained  by  solving  the  ray  equations  with 
initial  conditions 

S(*o)  -  ^vo'  ^(*o^  =  Pv  '  '^f^o^  =  ^<V^^  =  -•  ^o^V^^-  (^^5 

Here  P  =  (p^,...,p  )  is  arbitrary,  and  we  obtain  the  p-paramerter  family  of  jnye, 
Xy  «»  x^(tjP)  for  which  x  (t  |  P)  »  x  ,  (in  fact  we  obtain  two  such  families,  one 
corresponding  to  each  of  the  two  values  of  h).  When  the  rays  are  determined, 

s(t,X)  is  given,  on  each  ray  by 

t 
6(t)  -  8[t,X(tj  P)J  =  Sq  ■*■  /  >^[^C*'^  ^^*   ^^*'-  (^iJ 


Here  X  (X,K)  -  k^  hj^  -  h. 

In  solving  the  transport  equations  to  determine  the  functions  z^  "' 
along  the  rays  emanating  tram  M,  our  earlier  argument  is  valid  if  we  replace 
Z  «=  (0-,.,,,o  )  everywhere  by  P  =  (p-,,...,P  )  and  replace  t  in  (26)  and  (27) 
by  t,  ,  where  *-,  ^  *i  •  ^  course  now 
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a(x,...,x)  r^bc  (tjp)l 


functions  z^  '  become  singular  there,  which  is  what  we  expect  at  a  causti 


If  we  set 


^      z(\)\i(\)\'-'^ 


V*o      '-^hJ—  ^"^ 


then,  we  may  let  t,-t     In  (2  7)   (for  m  =  0).     Thus  we  obtain 

z(t)  =rc(t)|j(t)i-^/2.  (kh) 

Here  z(t)  =  z[t,X(t;P)J,  c(t)  =  c[x(t;P)],  j(t)  =  j(t,P),  w  =  h(X^,P),  and 
z"  =  z*  (P).  At  the  beginning  of  the  preceedlng  paragraph  we  stated  that  the  waves 
asSosciated  with  the  point  inarJ.fold  M  are  determined  by  the  values  of  s  and  z^^' 
on  M.  We  see  now  that  this  statement  must  be  modified,  for  the  z^''^'  are 


infinite  on  M.  In  fact,  the  leading  term  of  each  of  the  two  waves  emanating 
from 


from  M  is  uniquely  determined  by  the  values  of  the  number  s  and  the  function 


For  the  case  of  constant  coefficients,  k  is  constant  on  each  ray,  i.e. 
k^  =  k^(t^)  =  p^,  and  the  rays  are  given  by 

^v  =  ^o  -^  ^""-^K'   Sv  =  \  CP)  =  h^  (K)  (k3) 

V  V 

Purthennore  (kl)   becomes 
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s(t)  =  s^ 


+  (t-t^)i^    ;    ^^  =  k^h^.  h  =    /(K),  ik6) 


and 


j(t,K)  =  det   (t-t^)h^  ^       =   (t-^.)^  <i^(\  k   ^-  ('^^J 

'  V  i'  V  i 

Finally,  ve  may  replace   (kS)  by 

z(t)   =  c  zj  (t.t^)'^det  (h,^  3^    )  j-^/2.  (1.8) 

It  is  of  interest  to  note  that  the  discussion  of  this  subsection  is  not 
valid  for  a  non-dispersive  equation.    For  such  equations  the  rays  emanating 
from  a  point  all  lie  on  the  surface  of  the  characteristic  conoid  determined  by 
that  point,  hence  the  rays  do  not  fill  a  region  of  t,X- space  and  the  concept  of 
a  "tube"  becomes  meaningless.   Furthermore,  the  factor  det(h^  k. )  in  (k8)   vanishes 

in  the  non- dispersive  case. 

3  .  Initial-boundary  Value  Problems  for  the  Homogeneous  Equation 

In  this  section  we  will  show  how  siuns  of  one  or  more  waves  can  be  constructed 
to  yield  asymptotic  solutions  of  initial  and  initial-boundary  value  problems  for 
(1.17),  We  restrict  our  attention  to  the  simple  case  of  one  space  dimension  (|i=l) 
and  shall  obtain  just  the  leading  term  of  the  asymptotic  expansion.  Our  examples 

are  intended  to  be  only  illustrations  of  methods  which  are  applicable 
to  higher  dimensional  problems  as  well  as  to  other  differential  equations.   Therefore 
we  shall  avoid  using  special  methods  which  are  suitable  only  for  the  particular 
equation  under  consideration. 

We  consider  the  following  initial  value  problem; 


c^(x)u^.u^^-  aV(x)u  =  0,  (1) 

u(0,x)  =  u^(x;  -K),   Uj.(0,X)  =  Uj^(x;  7^).  (2) 


This  can  be  seen  by  differentiating  the  Euler  equation  k^  -^  =  h  with  respect 
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In  some  cases  we  shall  also  prescribe  a  homogeneous  boundary  condition.  An 
interesting  feature  of  the  asymptotic  method  is  that  entirely  different  types 
of  solutions  are  obtained,  depending  upon  how  the  parameter  A  appears  in  the 
initial  data  u , u .  We  consider  first  the  initial  value  problem  with  the  simplest 
type  of  initial  data. 
3.1  Oscillatory  initial  data 

Let  us  choose  initial  data  which  already  have  the  form  of  a  wave.  If 
v(t,x)  ^   z(t,x)e^^^^*^''^  then  ,^(t,x)  ~  iAs^ze^^^  =  -I7^(s^)2  e^^^   where 

h(k)  =  ^  h  (k)  =  tjc^^'  .  Let  us  therefore  consider  the  Initial 

value  problem  (l,2)  with 

uJx;A)  ^  z^(x)e^^o^^^  u^(x,A)  ^  .iAhrs;(x)Jz  Jx)e^^^o(^\  0) 

(Actually  we  are  simultaneously  considering  two  initial  value  problems,  since 
h(k)  has  two  values. )  The  functions  z  and  s  are  given. 

If  we  attempt  to  solve  this  inj.tial  value  problem  by  constructing  a 
wave, 

u(t,x)^z(t,x)e^^'^*'^^  {h) 

we  see  that  (3)  yields  the  values  of  s  and  z  at  t  =  0: 

s(0,  a)  =  s^(a)  and  z(0,a)  =  ^^i^)-  (5) 

Then  z(t,x)  and  s(t,x)  can  be  determined  by  the  method  of  section  2.  Let  us 
carry  out  the  details  for  the  case  of  constant  coefficients.  Prom  (2.I7) 


k^(a)  =  s;(a);  (.^(a)  =h[s;(a)].  (6) 
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We  then  solve  the  ray  equations  with  k  (a),  ^  (o)   given  by  (6),   t  (a)  =  0,  and 
X  (o)  =  a.  From  (2. 30)  and  (2. 31)  we  see  that  the  ray  emanating  from  the  iJoint 
(t,x)  =  (0,a)  is  given  by 

X  =  a  +  gt;  g  =  h'  [s;(a)]  ;  (?) 

while  on  this  ray 

k  =  s;(a),  0)  =  h|V(a)],  s  =  s^(a)  +  i  t;  X=  kg-w.  (8) 

Here  h»  =  dh/dk. 

To  determine  z(t,x)  we  first  con5>ute  the  jacobian, 

j(t,a)=g  =  i+t|f  =  ifth"(s;)s;;.  (9) 

Then  j(0,a)  =  1,  and  (2.35)  yields  the  value  of  z  on  the  ray: 

z(t)   =  zja)[l  +  th"(s;)s^]-^/2.  (10) 

Thus  we  see  that  the  asymptotic  solution  of  the  initial  value  problem  is  given 
parametricaHy  (with  parameter  a)  by 

u~z^(a)  |l  +th"ps;(a)Js^(a)  j  '^/\xp  jiA(^^(a)  +  |^;  h' (s;).h(s^  )jt)j 
X  =  a  +  h»  rSrt(cf)lt  ;  h(k)  =  +  h  (k)  =  -   ^ck+b^.  (ll) 


Noir 


hi (k)  =  c^(c\^+  b^)"^/\  J     h"(k)  =  cV(c\^+  b^)-3/2  >  0.  (12) 


Q.  .  ^  .  --  '  O' 

Let  us  assume  for  the  sake  of  definlteness  tliat 


sj  (a)  <  0.  (13) 
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Let  us  also  denote  the  solutions  corresponding  to  the  two  values 
h  =  -  h  by  u+  .  Then  for  u   ve  see  from  (l2)  and  (13)  that  the  jacobian 
(10)  vanishes  only  for  t<  0,  and  the  solution  formula  (ll)  for  u_^  does  not 
become  singular  for  t  >  0.  However  it  is  easy  to  see  that 


u  (t,x)  =  u^(-t,x) 


(Ih) 


hence  the  representation  (ll)  for  u  becomes  singular  for  certain  points  in  the 

half- space  t  >  0.  It  is  easy  to  show  that  these  points  lie  on  the  canst ic 

curve  or  envelope   of  the  ray  system.  The  singularity  is  a  defect  of  the 

asymptotic  solution  (ll)  and  must  be  repaired  by  considering  a  different  kind 

of  asymptotic  expansion  in  the  neighborhood  of  the  caustic.  The  rays  of  u 

are  illustrated  in  the  following  figure.  The  corresponding  picture  for  u_|_  is 

obtained  by  reflection  in  the  x  axis, 
t 


caustic 
It  is  interesting  to  note  that  if  m  =  dt/dx  =  g"^  is  the  slope  of  a  ray,  then 

from  (12)  |m|  >  c'  ,  i.e.  c  is  an  upper  bound  on  the  group  speed. 

Let  us  now  set 


z^(x)  ^ 


.  ^o^-) 


(15) 


^oW^^l 


u(t,x)  =  |[u_^(t,x)-ujt,x)]. 


(16) 
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Then  u  is  an  asymptotic  solution  of  (l)  and  we  see  from  (3)  that  it  satisfies 
the  initial  conditions 

u(0,x)  -0;   u^(0,x)  -  iVjx)e^^o(^\  (I7) 

Thus  the  asymptotic  solution  of  the  initial  value  problem  with  initial  conditions 
(17)  consists  of  a  sum  of  two  waves.     Finally,  let  us  set 

u(t,x)    =  U2Ct,x)   +  illt  ^1^*^^)  (1®) 


.  iAs.(X) 

u^(0,x)  -  0,  1^  u^(0,x)  ^    iAyj(x)e       ^       ',   3  =  1,2.  (19) 

Then  u  is  the  solution  of  (l)  which  satisfies  the  initial  conditions 

i7vs   (x)  iAs  (x) 

u(0,x)  ^  y3^(x)e       ^         ,   n^(0,x)  ^  ±-Ky^(x)e       ""       .  (20) 

(in  deriving  (20)  we  have  used  the  fact  that,  at  t  =  0,  u, .  .  =  c  u,   -  At)  u..  -^  0.  ) 
Thus  we  see  tlaat  the  asymptotic  solution  of  the  initial  value  problem  with 

oscillatory  initial  data  (20)  can,  in  general,  be  obtained  as  a  sum  of  waves. 
It  is  interesting  to  note  the  special  features  of  our  result  when  b  =  0 
In  this  case  the  parameter  A  is  absent  from  the  partial  differential  equation 
(1)  (in  fact  that  equation  becomes  merely  the  wave  equation)  and  appears  only 
in  the  initial  data.  Fran  (12)  we  see  that  the  group  speed  g  =  h' (k)  becomes 
g  =  +  c,  and  the  rays  lie  on  the  characteristic  surfaces  x  +  ct  =  constant. 
In  this  case  the  expsnion  reduces  to  Lax's  expansion  fh] . 
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3.2  the  phase-shift  at  a  caustic 

In  appendix  I  the  exact   solution  of  the  initial  value  problem  (l,2)    (with 
constant  coefficients)  is  given.      Let  us  set  h  =  +  h     in  (3)  and  examine  the 
solution  u     by  the  method  of   appendix  I.      From  (l.lO)  we  see  that 

where 

^^(k,o)  =  k(x-a);h^(k)t  +  sja).  (22) 

We  now  ajjply  the  stationary  phase  formula  (ll.  12)  of  appendix  II  to  obtain 
the  asymptotic  expansion  of  (21)  for  A  -^  «>.  The  stationary  point  is  determined 
by  the  conditions 

le  K  =  ^-^  '"■  K^^^^  =  °'  ^K  =  ^o^"^'^  =  °-  ^23) 

At  the  stationary  point  h  [s'(a)'|/h  (k)  =  1  therefore  the  term  corresponding  to 
the  lower  sign  in  (21)  vanishes.  At  the  stationary  point  (for  the  term  with  the 
upper  sign)  (22)  and  (23)  yield 

X  =  a  +  h;[s;(cT)Jt  ,  k  =  s;(a),  {2k) 
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ak       "^aa 


Hence 


det(^^.)  =.(1  +h;;  s;;t).  (27) 


It  follows  new  from  (II.12)  that 

"+-  z^<a)  ||1  f  th;rs;(a)]s;;(a)l]   -V2  exp|iA(s^(a)  +  [3ih^(-;)-hJs;)]t ) 

x=      a  +  h^[s^(a)Jt        ,     h^(k)  -/c^^-rt)^  .  (28) 

In  scsne  neighborhood  of  the  initial  plane  t  =  0,  1 1  +  th"  s"  |  =  l+th"s"  , 

and  we  see  that  (28)  agrees  exactly  with  (ll)  (for  h  =  +  h  )  provided  sig(<})  .)  =  0. 

In  order  to  determine  sig((}>  .)  we  conpjte  the  eigenvalues  r-,rp  of  (<}j  . ). 
They  are  the  roots  of  the  polynomial 

'^^^(^vj  -  ^^0^  =  (-h;;t.r)(sj.r).l  =  r^^  (h;t.s;)r  -  s^h^t-l.  (29) 

Thus 

2r^  =  (el  -  h;;t)+[(s^  +  h;;t)^  +  ^=^/^  (30) 

In  a  neighborhood  of  t  =  0  it  is  now  clear  that  r^  >  0,  r2<  0;  hence 

sig  (^  .)  =  sgn  r,  +  sgn  r^  =  0.  As  a  function  of  t,  sig(<|)  ,)  changes  discontinuously 

VJ  X  c.  Vj 

when  (27)  vanishes,  i.e.  when 
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t=t^(a)  =  .(h;;  s;;)-^  (31) 

We  recall  that  h"  >  0  and  s*  <  0.  Therefore  t(o)  >  0.  Prom  (30)  we  see  that  at 
00              o 

* '  V  (32) 


are  negative.  Thus 

3ig(Kj)  =  <;  ^.  (33) 

'^'^     1  -2     for   t  >  t^ 

The  result  (ll)  for  u  which  was  obtained  without  the  use  of  the  solution 
formula  (l.lO)  is  valid  in  the  region  t  <  t  (a).  We  now  see  that  that  result  is 
verified  by  (28)  and  (33).  However  (28)  and  (33)  are  also  valid  for  t  >  t  (a). 
This  suggests  a  rule  for  extending  our  general  method  to  the  region  t  >  t(2:), 
beyond  a  caustic  t  =  t(E):  The  formula  (2,28)  for  the  function  z  becones 
indeterminate  for  t  >  t(E)  because  j(t^)/j(t)  beccmes  negative.  The  indeterminacy 
is  removed  by  replacing  j/ j  by  |  j/ j  |  and  introducing  the  phase- shift  factor 
exp  \  -ijt/Sj-  .    This  "phase-shift  rule"  is  identical  to  an  analogous  rule  for  the 
reduced  wave  equation  [9]  . 

In  section  3.1  we  noted  that  the  asyn^rtotic  formula  f&ils  (i.e.  becomes  infinite) 
at  the  caustic  t  =  t(o).  The  correct  expansion  at  t  =  t(a)  could  be  obtained  by 
expejiding  the  integral  (21)  by  means  of  a  different  stationary  phase  formula  which 
is  valid  when  dert((})  )  vanishes  at  the  stationary-  point. 
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3.3  A  single  Initial-boundary  value  problem 

The  Initial  value  problem  for  u  ,  required  the  solution  of  (l)  with  the 
initial  conditions 

u(0,x)  ~  z  Jx)e^^o(^>,  ^  (0,x)  ^  i?yh^[sj(x)]z  Jx)e^'^o^^^  (3I;) 

The  asynptotic  solution  was  given  by  (ll)  with  h  =  -h  .  Let  us  now  modify  the 
problem  by  adding  the  boundary  condition  at  x  =  0, 

u(t,0)  =  0.  (35) 

The  solution  in  the  quarter-space  t  >  0,  x  >  0  is  required.  For  consistency 
ve  must  assinne  that  z  (o)  =  0. 

The  asymptotic  solution  of  this  initial-boundary  value  problem  is 
obtained  by  adding  to  u  a  reflected  wave  u  (t,x)  which,  together  with  u 
satisfies  (35).  The  reflected  wave  must  vanish  at  t  =  0  so  that  (3^+)  is 
unchanged,  Tbaxe  we  set 

u  =  u_,  +  u^,  u_(t,x)  ~  z(t,x)e^'^^*'^^  u^(t,x)  -^   z^(t,x)e^'^r^*'^^       (36) 

Here  z  and  s  are  given  along  the  "incident"  rays  (?)  by  (lO)  and  (8). 
By  Inserting  (36)  in  (35)  we  obtain 

8^(t,0)  =  s(t,0),  (37) 

Zr^*'°^  =  -z(t,0),  (38) 


-3^- 


and  differentiation  of  (37)  with  respect  to  t  yields 

a)^(t,0)=  u)(t,0).  (39) 

2    2. 2  2 
Since  "both  incident  and  reflected  waves  satisfy  the  dispersion  equation  w  =  c  k  -fc 

we  conclude  that  k  (t,0)  =  -k(t,0),  hence  g  (t,0)  =  -  g(t,0).  In  order  that  the 

reflected  rays  should  not  intersect  the  initial  line,^  t  =  0,  x  >  0,  we  must  take 

k^(t,0)  =  -k(t,0),  g^(t,0)  =  .g(t,0).  (ko) 

It  is  convenient  to  express  the  boundary  pai-ametrically  in  terms  of 
the  same  parameter,  a,   used  in  describing  the  incident  wave  u  ,  From  (7) 
we  see  that  at  x  =  0,  a  +  g(a)t  =  0.  Hence  the  boundary  is  given  parametrically 
by 

X  =  x^(a)  =0   ,  t  =  t^(a)  =  -a/g(a).  (kl) 

Then,  on  the  reflected  ray  emanating  from  the  boundary  at  [t  (a),Oj,  (39)*  (^)t 
and  (6-8)  yield 


=  -^oC^J^'^JJ'  \  =  -k  =  -sj(a),  g^  =  -g  «  hj[s'(a)],  i^  =  k^g^.  u^ 
=  kg-u  =  (/ 


(^) 


From.   (2.31)  the  parametric  equation  of  the  reflected  ray  is 

X  =  gj.(t-t^)  =  -g(t-t^)  =  -a-gt  =  -o  +  th«[s^(a)J,    (k3) 

while  (2.32),  (37),  (8),  and  (k2)   yield 

^  =  s[t^(a)]  +  ^j.(t-t^)  =  s^(0)  +  J^t^   +  X(t-t^)  =  s^(ct)  +  2t,  (kk) 


*This  condition  is  clearly  required  if  u  =  u  +  u  is  to  satisfy  the  initial 
conditions,  for  then  we  must  require  that  u  =  0  at  t  =  0,  x  >  0. 


-35- 


From  (U3)  we  may  now  compute  the  jacdblan. 


|^i.th;[sj(c)]s;;(a)}    . 


Then  (38),  (2.28),  and  (lO)  yield 
1/2 


z^(t)=z[t  (a)J 


.z^(a)    [l-t^h;  sj]-V2 

.zja)[i-t  hj  b;;]-V2. 


1-t X  s"     ,1/2 

00      O      v»  ' 

l.th  "s" 


(h5) 


(1*6) 


and 


Summarizing  our  results,  we  see  that  u  *  u     +  u    where 
u^».  -z^(0)[i-th^  ^o]~^^^    ^"^  {  ^'^^%  +  -^  t)  j  ,  X  -  -a-gt 


u_^  z^(a)[l-th^  3^""^^^  ^^  f  ^^^%  +it)  I    ,  X  =  a-tgt. 


(^7) 


(^) 


This  con5)letes  the  solutlcm  of  the  initial-boundary  value  problem. 
From  (47)  and  (1*8)  we  note  that 


u^(t,x)  =  -u_(t,-x). 


(k9) 


Of  course  {k9)   is  easily  obtained  by  the  "image  method**  which  yields  the  solution 

of  the  InitiaX-boundary  Value  problem  at  once  in  terms  of  the  solirtion  u  of  the 

initial  value  problem.  Thus  this  method  provides  a  check  on  the  construction 

we  have  vised. 

The  above  problem,  though  essentially  trivial  suffices  to 

*  In  this  case  the  image  method  would  consist  in  extending  tlB  initial  data  to 
X  <  0  in  such  a  way  that  u(0,x)  and  u.  (0,x)  are  odd  function  of  x. 
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illustrate  the  general  procedure  for  the  asymFrtatic  solution  of  initial-boundary 
value  problems:  We  first  obtain  the  asymptotic  solution  of  the  initial  value 
problem  in  terms  of  one  or  more  "aves  emanating  from  a  space- like  initial  surface. 
Vhen  these  waves  are  "incident"  on  a  time- like  bcSundary,  additional  reflected 
vaves  are  produced.  Rich  reflected  ^-ave  is  uniquely  determined  by  the  incident 
wave  and  the  boundary  condition.   Of  course  we  must  expect  complications  to 
occtxr  when  incident  rays  are  tangent  to  the  boundary  or  when  boundaries  are  not 
smooth.   In  such  cases  there  will  be  "shadows"  and  we  shall  have  to  study 
"diffracted"  waves  In  the  shadow  region.* 

The  incident  and  reflected  rays  of  our  simple  problem  are  illustrated  in 
the  follor./ing  figure: 


3.^  rapidly  varying  initial  data 

Our  results  so  far  have  been  restricted  to  problems  with  oscillatory  initial 
data.  In  this,  and  succeeding  sections  ve  shall  consider  other  kinds  of  initial 
data.   In  each  case  the  initial  data  give  rise  to  one  or  more  ../aves.   Then  if 
the  problem  involves  boundaries,  reflected  --aves  will  be  produced  just  as  in 
section  3.3.   Thus  the  essential  new  features  are  embodied  in  problems  "ithout 


*Such  diffracted  waves  for  elliptic  problems  are  exhaustively  treated  in  [ij  and 
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boundaries  and  we  may  restrict  our  attention  to  pure  initial  value  probelms. 

In  this  section  we  consider  the  initial  value  probelm  for  equation  (l) 
with  initial  data  of  the  form 

u(0,x)  =  u^(x|  A)  =  u^(?oc),  u^(0,x)  =  u^(x;7v)  =  7vu^(>x).  (50) 

Such  dataare  called  rapidly  varytng  because  of  their  large  derivatives.  Inoterms 
of  the  transformation  (l.  3 1)  th->  s  problem  is  equivalent  to  the  problem  for  the 
differential  equation  without  a  parameter, 

(=')^«'x.x.-'''t.f  -("''""'  =°  (51) 

and  with  initial  data  of  the  form 

u«(0,x')  =  u^(x'),  u'^.(0,x»)  =  u^(x»),  (52) 

i.e.  also  independent  of  the  pai-ameter.  Thus  we  see  that  rapidly  varjring  initial 
data  are  quite  natural  for  certain  problems. 

We  assume  that  the  functions  u  (x*),  u.(x*)  have  cocyact  support,  i.e. 
vanish  outside  a  finite  region  called  the  support  of  these  functions.  Then  for 
X  -♦•  the  support  of  u  (>«)  and  Ao,  (?oc)  shrinks  to  the  point  x  =  0.  This  suggests 
that  the  waves  produced  by  the  rapidly  varying  initial  data  (50)  emanate  from  the 
single  point  (0,0).  In  section  (2.7)  we  saw  that  a  point  initial  manifold,  in 
general,  gives  rise  to  two  waves.  Therefore  we  set 

iTvs 


u  »  u  +u  ;  u.  -  z   e   +  ,  (53) 
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The  functions  z   and  s  are  given  by  the  method  of  section  2.7  except  that  the 

numbers  s   and  the  functions  z*  (p)  are  undetermined.  In  order  to  determine 

the  rave  proceed  as  follows:  Plrst  we  complete  the  construction  for  the  case  of 

constant  coefficients.  Then  we  shall  compare  the  resiilt  with  the  asymptotic 

expansion  of  the  exact  solution  <>f  our  problem.  This  will  determine  s    and 

Z' .(p).  We  then  postulate  that  the  newly  determined  quantities  have  the  same  values 

for  the  case  of  variable  coefficients*. 
From  (2.U7)  the  rays  are  given  by 

X  =  gt;  g  =  h'(k)  =  -  h^(k).  (^h) 

Then  from  (2,1+8)  and  (2.50)  we  see  that  s  and  z_^  are  given  along  the  rays  by 

s+Ct)  =  %+-^  i^}     ^=  kh'(k)-h,  (55) 


and 


.^(t)  ='^^(k)c|t  h;;(k)ri/2  (56) 


Therefore  for  t  >  03u(t,x)  is  given  isarametrically  in  terms  of  the  i^rameter  k  by 

u^  ^      cz_^(k)lth;(k)|-^/2  exp  ^  i-K[s^^  +  (khj  -  h^)tj]  i 

(57) 

x  =  +  h^(k)t. 
The  exact  solution  of  the  initial  value  problem  (l,50)  with  constant 

coefficients  is  f^iven  by  (I.3)  with  a  given  by  (I.9).  We  therefore  insert  (50) 


*This  is  the  method  of  comparison  with  the  asymptotic  expansion  of  the  solution 
of  a  "cannonical  problem".   The  validity  of  the  procedure  for  the  problems 
treated  here  has  been  confirmed,  as  pointed  out  in  the  introduction. 
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in  (1.9).     This  yields 

(58) 


We  note  that  a     is  independent  of  7v.     Frcm  (l.3)^ 

u(t,x)  =  y~   y  a^(k)e^^^^W;       <|)  =  kx-ht;  ^  =  *o  =  -[^^"^  -  (59) 

+ 

Me   now  apply  the  one-dimensional  stationary  phase  formula  (II.I5)  to  each  of  the 
two  terms  in  (59).  The  stationary  point  is  determined  by 

^'(:0  =  x-gt  =  0,  g  =  h»(k)  =  -  hj^(k),  (60) 

and  at  the  stationary  point 

^•(k)  =  -h"(k)t  =  4h;;(k)ti  <})  =  (kg-h)t  =  +  (kh«  -  h^)t.  (61) 

Therefore  (11.12)  yields,  for  t  >  0 

^-  H     H^^^^  %(^^^      ^  |th;;(k)r^/^  exp|iAp(khJ  -  h^)t]  j  ; 
+  " 

X  =  ^   h»(k)t. 


(62) 


(Here  we  have  used  the  fact  that,  h"(k)>  0,  hence  sig  (()"  =  sgn  ^"   =  +1  for  t  >  0. ) 
We  now  note  tlmt  (57)  and  (62)  agree  exactly,  provided  we  set 

So+  '  0}       z^(k)  =  (-^)  /  —3 e      .  (63) 

Here  a^(k)  are  the  fourier  transforms  of  the  initial  data,  given  by  (58).  Thus 
we  have  verified  our  general  method,  and  have  determined  s  +  and  z^   . 
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The  asymptotic  solution  of  the  initial  value  jiroblein  for  the  case  of  variable 

coefficients  c(x),  "b(x)  can  novr  be  comp3.eted  by  the  method  of  section  (2.7)  with 

s  and  z'_j_  (p)  given  by  (63).   Of  course,  it  is  clear  that  in  that  case  we  should 

replace  c  and  b  in  (63)  by  c(0)  and  b(0).*  The  details  of  the  construction  will 

be  omitted. 

For  the  case  of  constant  coefficients,  the  rays  assosciated  with  each 

wave  are  illustrated  in  the  following  figure. 

t 


x=-ct 


>  X 


From  (12)  we  see  that  for  -00  <  k  <  »,  _c  <  h*(k)  <  c  hence  the  rays  all  lie  within 

the  characteristic  cone  x  =  -  ct  emanating  from  the  origin.  The  origin  itself 

is  a  caustic  of  the  ray  family  and  the  asymptotic  solution  is  singular  there. 

From  (12)  we  note  that  [h^  (k)]""^'^  -0(|kp/^)  as  k  -»  -  00,  therefore  the 

asymptotic  solution  also  becomes  singular  on  the  characteristic  cone  unless 

a,  (k)  =  0(|k|"-^'  ).   Singularities  of  asymptotic  solutions  on  caustics  and  other 

special  manifolds  (such  as  shadow  boundaries)  are  well  known  phenomena  in  the 

partial 
asymptotic  theory  of  differential  equations.  The  singularity  on  the  characteristic 

cone  appears  to  be  another  phenomenon  of  this  type. 

For  the  case  of  variable  coefficients  the  ray  picture  will  oi:'  coiirse  be 

similar  to  the  figure  except  that  the  rays  will,  in  general,  be  curves  lying 

within  the  characteristic  conoid. 


*(63)  depends  on  b  as  well  as  c  because  h  (k)  = /c  k"-ft»   appears  in  (58). 
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As  a  special  case  of  otir  method  we  obtain  the  asymptotic  expansion  of 
the  Rlemann  function  for  (l)  which  is  a  solution  of  the  initial  value  problem 
with 

u(0,x)  =  0   ,  u^(0,x)  =  S(x).  (6U) 


Since  5(x)  =  X6(Xx),  we  may  take  u  (x' )  =0  and  u^(x' )  =  6(x')  in  (50).  Then 


(58)  yields 


a^  (k)  =:^i[U«hJk)j-^  (65) 


and  (63)  becomes 

s^  =  0,      z^(k)   =  [I(2rt  X)V2  ^  hjk5]-l  e-  ^  '^/^  (66) 

3.5  Slowly  varying  initial  data 

In  view  of  our  definetion  (50)  of  "rapidly  varying"  Initial  data, 
it  is  appropriate  to  refer  to  initial  data  of  the  form 

u(0,x)  =  u^(x;  X)  =  u^(x),  u^(0,x)  =  u^(x;  X)  =  Xu^(x)         (6?) 

as  slowly  varying.  In  terms  of  the  transformation  (I.31)  the  initial  value 
problem  (1,67)  is  equivalent  to  the  initial  value  problem  for  (51)  with  initial 
data 

u'(0,x')  =u^(xVx),  u'^'(0,x')  =  u^(x7x).         (68) 
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The  asymptotic  solution  of  the  problem  {l,6l)   is  easily  obtained  by  observing 
that  (67)  is  a  special  case  of  oscillatory  initial  data(20),  with 
s,  =  Sp  =  0,  y,  =  u  ,  and  iyp  =  u, .  As  we  saw  in  section  (3.I)  the  asymptotic 
solution  of  such  problems  consists  of  a  sum  of  waves,  each  of  which  is  generated 
by  data  of  the  form  (3).  If  we  set  s^(x)  =  0  in  (3)  it  becomes 

u^(x;X)^  z^(x),   u^(x,X)  -V.  -iXh(0)E^(x);  h(0)  =  -  |b(x)|.    (69) 

For  the  case  of  constant  coefficients,  the  asymptotic  solution  can  be  obtained 
from  (11)  by  setting  s  =  0.  This  yields 

u^  (t,x)  ^  z^(x)exp{+  iX|bIt).  (70) 

The  rays,  in  this  case  are  straight  lines  parallel  to  the  t-axis,  i.e.  the 

group-speed  is  zero.  If  we  follow  the  procedure  outlined  in  (15-20)  we  obtain 

the  asymptotic  solution  of  the  initial  value  problem  (1,67)  for  constant  coefficients; 

u,(x) 
u(t,x)  ^u^(x)  cos(X|b|t)  +  ^U-  sin(X|b|t).  (7I) 


o 
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3.6  disconrtinuous  initial  data 

It  is  a  characteristic  featvire  of  as3niiptatic  solutions  of  differential 
equations  that  discontinuities  of  the  data  give  rise  to  additional  terms  in  the 
asymptotic  expansion.  In  order  to  illustrate  this  phenomena,  let  us  suppose  that 
the  function  z  (x)  in  the  initial  value  problem  (1,2,3)  is  sm.ooth,  except  for  a 
jump  discontinuity. 
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at  X  =  X  .  Pbr  definiteness  we  shall  take  h  =  -Hi  in  (3). 
o  o 

We  ass\me  that,  in  addition  to  the  wave  (ll)  generated  by  the  initial  data, 

two  secondary  waves  u  /s>  z_j_  e   +  emanate  from  the  point  of  discontinuity  (0,x  ) 

on  the  initial  line  t  «=  0   The  functions  z  and  s   are  given  by  the  method  of 

section  2.7  except  that  the  numbers  s  +  and  the  functions  z'.(p)  are  undetermined. 

o-  + 

In  order  to  determine  them  we  proceed  as  in  section  3.4.  The  first  step  is  to 
complete  the  construction  for  the  case  of  constant  coefficients. 
From  (2,47)  the  rays  are  given  by 

X  =  x^  +  gt;   g  =  h'(k)  =  \M.  (73) 

Prom  (2.U8)  and  (2.50),  s_^  and  z^  are  given  along  the  rays  by 

s+  (t)  -  s^^  +  it;  J(=  kh'-h  =  +  (kh^  -  h^)         (Ik) 


and 


Thus 


z+(t)  -r.^(k)c|th;;  (k)ri/2.  (75) 

u^-  cz^(k)|h;;(k)t|-^/^  expfiXJ-s^^  _+(kh;  -  h^)t]|  ; 


X  =  X  +  h'(k)t. 
o-  o^ 
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solution,  which  is  given  by  (21).   The  leading  term  of  the  asymptotic  expansion 
of  the  double  integral  (21)  comes  from  the  neighborhood  of  the  stationary  point 
and  has  already  been  computed  (see  (28,  33)-)   It  can  be  shown  that  the  next 
term  in  the  asymptotic  expansion  comes  from  the  neighborhood  of  the  discontinuity 
line  a  =  X     of  the  integrand  and  can  be  determined  as  follows: 

We  first  divide  the  integral  (21)  into  two  parts,  w^  corresponding  to 
the  strip  x^-€  <  ^  <  ^^  *  ^)   a"d  v-^   corresponding  to  the  complementary  region 
of  the  o,   k- plane.  We  choose  €  small  enough  so  that  the  stationary  point  lies  in 
the  domain  of  w  .  w  may  be  written  in  the  form 


where 


and 


r^     =  f     v^(k)e^'E-^o(^^)*l 


(77) 


dk.  (78) 


X    r  °  ixrs_^(o)-k?] 


da; 


q  ±(a,k)  -   ^  -    h  (k) 


(79) 


Next  we  divide  each  integral  (79)  into  two  parts  corresponding  to  the  intervals 
X  -  e  <  (J  <  X  2nd  X  <  a  <  X  +  e,  and  integrate  each  by  parts.  It  is  easily 
shown  (by  a  second  integration  by  parts)  that  the  resulting  integrals  are 
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of  order  l/"^   while  the  boundary  tennsare  of  order  1.   Of  the  four  boundary  terms, 

the  two  which  come  from  the  points  x  -  e  and  x  +  €  must  be  ignored  for  they 

are  exactly  cancelled  by  boundary  contributions  from  v    Then  the  boundary  terms 

he  points  x  -  0  and  x  +0  yieL 
o        o 

-fzl  q+(x  ,k)   iXfs  (x  )-kx1 
We  now  insert  (80)  in  (78),  This  fields 


where 

<f>^(k)  =k(x-x^)±h^(k)t.  (82) 

Finally  we  apply  the  stationary  phase  formula  to  (81 ).  The  stationary  point  is 
determined  by 

^»(k)  =x-x^±h;(k)t  =0,  (83) 

and  at  that  point, 

f  =  -+h;;(k)t,     ^  =  ±jkh;  -  h^t.  (ek) 

Then,  for  t  >  0,  (II.12)  yields  w^^v^u^  where 
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^         i[z Je^  "/''  (l±h„[s'(x^)]/h^(k)]       ,/.,     ,   p  -|, 


(85) 


X  -  x^  ±  h;(k)t. 


Since  the  integral  v^  contains  the  stationary  point,  its  asynrptotic 
expansion  is  again  given  by  (28,33),  and  coiresponds  to  the  primary  vave.   Therefore 
the  two  secondary  vaves  correspond  to  (85).  If  we  nov;  compare  (65)  and  (76)  we 
see  that  they  agree  exactly  provided  we  set 


V>-o(V'  ^'.(■''"i^^i^T^^'c 


i[ge«"/'*     (i±hjs;(x„)]/hjic)} 


6'(x  )-k 

O^  O 


(86) 


VJe  note  that  the  primary  wave  (28)  is  of  order  1  while  the  secondary 

-I/2 
waves  are  of  order  X  '  ,  The  rays  of  the  prinary  wave  are  illustrated  in  the 

figure  of  section  3.1.  Those  of  the  secondary  waves  are  illustrated  in  the 

t 

following  figure: 


■Caustic  of  primary  -v-ave 
They  all  lie  within  the  characteristic  cone  (dashed  lines)  emanating  from  the 

point  (0,x  ).  One  ray  (dotted  line)  coincides  with  a  ray  of  the  primary  wavie. 

Prom  (7)  we  see  that  the  eqution  of  this  ray  is 


^ "  ^o  -^  h;r3;(x^)jt. 


(87) 


-1^7- 


Prom  (73)  we  see  that  the  ray 


X  =  x^  +  h^(k)t  (88) 


of  the  secondary  wave  u   coincides  with  it,  provided 


.;Cx^).  (89) 


But  for  this  value  of  k  we  see  from  (86)  that  z'.(k)  is  infinite   (while  z'    is 

indeterminate).     Thus  the  secondary  wave  u^    is  infinite  along  the  ray  which 

~l/2  ^/'' 

co-lncides  with  a  ray  of  the  primary  wave.      Since   fh"   (k)]  '      =  0(|k|   '■")     as 

k  -» 00,  it  is  easily  seen  that  both  secondary  waves  become  singular  on  the 

characteristic  cone. 

As  in  section   (3.'+),  the  formulas   (86)   for  s  ±  and  z'_^(p)  can  be  used 

to  complete  the  construction  of  the  secondarj'-  waves  u^     for  the  case  of  variable 

coefficients.      In  doing  so  we  must  take  c   =  c(x  )  and  b  =  b(x  )  in   (86).      Once 

u^  are  determined,  the  asymptotic   solution  of  the  initial  value  problem  (l,2,3) 

(with  z     discontinuous  at  x  =  x      and  h  =  +h   )  is  given  by 

u  =  u^  +  u^  +  u_    +  0(l/X).  (90) 

Here  u  is  given  by  (ll)  for  the  case  of  constant  coefficients.   For  the  case 

of  variable  coefficients  u  can  be  determined  by  the  method  of  section  2,   using  (86). 
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k.      Initial-boundary  value  problems  for  the  Inhomogeneous  etjuation 

In  this  section  we  shall  develop  methods  for  obtaining  the  asymptotic 
solution  of  initial-boundary  value  problems  for  the  inhomogeneous  equation 

c-(x)u^-u^^-  xV(x)u  =  x2f(t,x;X).  (1) 

VJithout  loss  of  generality,  we  may  assume  that  the  initial  data,  which  we 
prescribe  at  t  =  t,<  0,  is  homogeneous,  i.e. 

u(t^,x)  =  u^(t^,x)  =  0.  (2) 

(The  solution  of  the  problem  for  (l)  with  inhomogeneous  initial  data  is  then  given  by 
u,  +  Up  where  u.  is  the  solution  of  (1,2)  and  u^  is  a  solution  of  the  problem 

(3.1,2)  for  the  homogeneous  equation.). 

2 
We  shall  assume  that  the  source  term  X  f  gives  rise  to  one  or  more 

waves.  The  sum  of  these  waves  provides  the  asymptotic  solution  of  the  initial 

value  problem.  For  initial-boundary  value  problems  we  assume  that  a  secondary 

reflected  wave  is  produced  (as  in  section  3.3)  when  a  primary  wave  is  incident  on 

the  boundary.   Since  we  have  already  examined  this  phenomenon, at  least  in  a 

simple  case,  we  shall  restrict  our  attention  here  to  problems  without  boundaries. 

In  the  following  subsections  we  shall  construct  the  waves  produced  by 
several  different  types  of  sources. 
k.l     rapidly  varying  source 

The  sovirce- function  will  be  said  to  be  rapidly- varying  if  it  is  of  the 
form 

X^f(t,x;X)  =  x2g(Xt,Xx).  (3) 
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In  terms  of  the  transformation  (l.lo)  the  initial- value  problem  (1,2,3)  is  then 
equivalent  to  the  problem 

(c')2u'^,^..u'^,^,  -  (b')2u'  =g(t',x'),  (U) 

u'(U^^x')  =u'^,(H^^x')  =0.  (5) 

Hence  we  see  that  this  type  of  source  is  quite  natural  for  the  problem  without 
a  parameter  in  the  differential  equation  and  with  homogeneous  initial  conditions 
in  the  "distant  past",  i.e.  at  t'=  H.,  (t.<  0,  X  -»c»). 

Let  us  assume  that  the  function  g(t',x')  has  compact  support.*  Then 
for  X  ^co  the  support  of  g(Xt,Xx)  shrinks  to  the  point  (t,x)  =  (O^O).   Therefore 
ve  assume  that  the  source  (3)  gives  rise  to  craves  emanating  from  this  point.   If 

we  set 

iXs^ 
u+  ~  z^e    ,  (6) 

then  the  functions  z+  and  s^  are  given  by  the  method  of  section  (2.7)  except  that 
the  numbers  s  +  and  the  functions  z^(p)  are  undetermined.     As  in  section  (3.''+), 
we  shall  carry  out  the  construction  for  the  case  of  constant  coefficients,  and 
then  shall  determine  s    ,  and  'z_^_  by  comparison  with  the  asymptotic  expansion  of  the 
exact  solution. 

The  asymptotic   solution  for  the  case  of  constant  coefficients  is  the 
same  as  that  of  section   (3.^)  exceprfc  for  the  values  of  s   ^  and  z_^_.      Tlierefore 
u+  is  given  along  the  rays 

x  =  ±  h;(k)t  (7) 

"♦  Actually  it  is  sufficient  to  assume  that  the  decay  at  infinity  of  g(t',x'")  is 
sufficiently  rapid  for  the  Integral  (11 )  to  exist. 
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by 

u^  »  cz^(k)|th;;(k)|-^/2  e::p  [  iX[s^^±(kh;  -  h^)t]j  -  (8) 

The  exact  solution  of  the  initial  value  problem  (1,2)  is  given  by  (I. 23) 
of  appendix  I.  We  insert  (3)  in  that  equation  and  then  introduce  the  transformation 
x'  =  Xy,  t'  =  Xs.   The  result  is  /_» 

00      00        \t 

/+   ^^   T —   1   r  ^1  r  ^  1  r  ^4.1  gCtSx')  -i[kx'-ht'1-HXrkx-htl. 
u(t,x;X)=^  ^   J^^^J^^     dx'^^^  dt'S^^^e  L      J   L    J 

Here  h  =  ih   (k).      For  X  -» «  it  is  clear  from  (9)  that 


"  ~  I!    ^   I    ^    e^HO')  aK,  (10) 

00  00 

a(k,h)   =   r   dx'    Tdt'   g(f,x')e-^P'^'-^*'^,  (11) 

^00  =  kx-ht.  (12) 


where 


and 


We  now  apply  the  method  of  stationary  phase  (appendix  II)  to  (10). 
The  stationary  point  is  determined  by 

<})'(k)  =  x-h't  =  0,  (13) 

and  at  the  stationary  point, 

^  ==  (kh'-h)t  =  +(kh^  -  h^)t;  f  =  -h"t  =  '^'^  t.  (lU) 


♦Actually  it  is  sufficient  to  assvune  that  the  decay  at  infinity  of  g(t;  x')  Is 
sufficiently  rapid  for  the  integral  (11 )  to  exist. 
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According  to  (3.12)  h"  >  0  therefore  f or  t  >  0  sgn  ^"   =    '^  1.     It  follows  from 
(11.12)  that  u  .X.  u^  +  u  ,  where 

jj^         |th;;(k)|-^/-  exp|^iX[^±(kh;  .  h^)tj+  31   n/u]  ; 


2(2nX)"/"h^ 


(15) 


X  =  ±  h^(k)t. 


This  eqimtion  agrees  exactly  with  (8)  if  we  set 


,  =  0,  ..(.)  =  Ihtlo^   ,«i./S  ,  =  ^VT^a 


2(2nX)'^ch^(k) 


(16) 


Here  a(k,h)  is  given  by  (11 ).   (l6)  can  be  used  to  complete  the  asymptotic 
solution  of  the  problem  with  variable  coefficient  functions.  In  so  doing  we 
must  set  c  =  c(0)  and  b  =  b(0)  in  (l6). 

Two  special  cases  of  the  results  of  this  section  are  of  some  interest. 
One  is  the  point  source  with 

g(t' ,x')  =g(f)5(x').  (17) 

The  other  is  the  source 

g(tSx')  =  S(f)5(x')  (18) 

which  determines  the  Green's  function  of  the  initial  value  problem. 
For  the  point  source,  we  see  from  (ll)  that 

00 

a(k,h)  =a(h)  =  f   g(t')e^^*'  df,  (19) 

while,  for  the  Green's  function 
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t(k,h)  »     1.  (20) 


i<-.2     Source  with  rapid  time  variation 

Let  us  now  consider  a  source  function  of  the  form 

f  (t,x;X)   =  g  ^  >^[t-t^(x)J  ,x  }    .  (21) 

If  the  function  g(t',x)  has  compact  suprport*  then  as  X  -»(»  the  supi^rt  of  f 
shrinks  to  the  cxirve  t  =  t   (x).     In  order  to  conform  to  the  notation  of  section  2, 
we  express  this  curve  in  the  parametric   form 

t  «  t^(a)    ,     X  =  x^(a)   5  a.  (22) 

To  simplify  our  subsequent  work  it  is  convenient  to  assume  that  t  (a)  is  a  monotonic 
function,  and  for  definiteness  we  shall  take 

t^C")  >  0.  (23) 

VJe  assume  that  two  waves  are  produced  by  the  source   (21)  and  that  they 

emanate  from  the  curve   (22).     V7e  also  assume  that  on  this  curve,  s(t,x)  is  constant 
for  each  wave,   i.e. 

s[t^(a),a]   =  s^^  =  const.  (2U) 

Differentiation  of  (2l|)  with  respect  to  o  yields  the  condition 

-cjt.^  +  k  =  0  or  w/k  =  1/t^  «  ~.  (25) 

(The  condition  (25)  may  he  stated  as  follows:     The  "phase  speed"  w/k  is  equal  to 
the  speed  dx/dt  of  the  source.).     Bj'  squaring  (25)  we  obtain  k     =  (t')  h  (k,a), 

*  It  is  sufficient  to  assume  that  the  integiral  {hk)  exists. 
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2 
and  since  h  depends  on  k  only  as  a  function  of  k  this  equation  has  two  solutions 

k  =  ^Jo).  (26) 

(By  using  the  formula  (2£)  for  the  hamiltonian  h  it  is  easy  to  find  the  solutions 
(26)  explicitly.   However,  since  we  wish  to  keep  our  methods  as  general  as  possible 
we  shall  not  make  use  of  the  explicit  expression  for  k  . ). 

We  now  impose  the  condition  that  the  solutions  (26)  be  real,  and  we  take 
k  (o)  >  0.   For  the  case  of  equation  (l)  it  is  easy  to  show  that  this  condition 


is  equivalent  to 


|c(a)t;(a)|  <1   or  \c\   <  |g|,  (27) 


i.e.   the  source  speed  must      be  greater  than  the  characteristic  speed  Ic(, 
(Geometrically,  the  condition  is  that  the  curve  (22)  should  be  space- like. ). 
From  (25)  and  (2.6) 

u)  =  ioja)  =  ik^(o)/t^  =  ±h^[k^(0),a]  (28) 

and  from  (2.8) 

bh  Sh 

The  value  of  the  function  z  for  each  wave  is  determined  by  its  value 

zrto(o),cT]  =  z^a)  (30) 

on  the  curve  (22).  If  z_^{a)   and  the  constants  s  ^  were  known,  the  construction 
of  the  two  waves  could  be  completed  by  the  methods  of  section  2.  We  shall 
determine  z^  and  s^^  by  completing  the  construction  for  the  case  of  constant 
coefficients  and  comparing  the  result  with  the  asymptotic  expansion  of  the  exact 
solution. 
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If  the  coefficients  are  constant,  the  rays  are  given  by 

X  =  a  +  g^(a)  [t.t^(a)]  =  a   +  h;[k^(a)J  [t-t^(a)]  (3I) 

and  s  is  given  along  the  rays  by 

s  =  Sq+  +  (t-t^),  (32) 

where 

;=ks.c.  =  ±[kja)h;(k^)-h^(k^)].  (33) 

Erora  (2,35)  and  (30)^  z  is  given  along  the  rays  by 

^(t)   «  z^  (o)  j](t^)/j(t)]V2^  (3i^j 

where  j(t,cr)  =  dx/dcr.  In  order  to  calculate  this  derivative,  using  (3I),  we  shall 
need  the  value  of  k'(cf)  =  dk°/do.  This  is  obtained  by  differentiating  icplicitly 
the  equation 

\KK   =  K  (35) 

which  is  a  consequence  of  (25),  (26),   and  (28).  The  result  is 

^;(^)=Vo"/(i-h;^;)-  (36) 


Then 


j(t,a)  =  dx/da  =  i-h^t;  +  (t-t^)h;;  h^tyd-h^t;).  (37) 


It   follows  that 

j(t^,a)=i.h;t;,  (38) 

*  We  note  that  6A^)  is  the  same  for  both  i/aves,  therefore  the  rays  are  the  same. 


-55- 


and  (3k)   yields 

^(t)  =  z^(,)  [(i-h;t;)^  -H  (t-t^)h;;  h^t;]-V2  (39) 

If  we  denote  the  two  waves  by  u,  and  u  ,  then  the  asymptotic  solution  of  the 
initial  value  problem  (1,2,21)  is  given  by  u  =  u  +  u  ,  where 

u^  ^  ^±(-)[a-h;t;)^  Mt.t^)h;;h^t^-V2  exp[ix[s^,±(k^h;.h^)(t.t^i]j  ; 
X  =  a  +h^(t-t^).  (^0) 

Here,  the  argument  of  h  ,  h',  h"  is  k  (a). 

The  exact  solution  of  the  problem  is  given  by  (I.23),  We  insert  (21) 
in  that  equation,  change  y  to  cr,  and  introduce  the  transformation  t'  =  X^s-t  (cr)^ 
in  the  inner  integral.  This  yields 

00  00  X   ft-t   (a)1  (Ul) 

L  O  J  ^/.,       _N  r  ^-  -I 


X 

lijti 


/     dJ:  /       da     /  dt'       S(^';^)     exp  i'lxrk(x-(7)-h(t-t    )"|-»lhf 


Here  h  »  ih  (k).  We  assuiae  that 


V  ^o^''^^  *•  ^^^^ 


Then,  for  X  -» », 


where 


Z_   '^"^  ^.00      ^.co       h(k) 


t(h,a)  =  f     g(f,ci)e^^*'df  (hk) 
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<t)(k,a)  =k(x-a)  .  h[t-t^(a)J.  (1^5) 

VJe  now  apply  the  method  of  stationary  phase  to  (U3).     The  stationary 
point  is  determined  by  the  conditions 

<j>j^  =  x-a.h'[t-t^(a)]    =0,       ^o  =  -k  +  ht;(a)   =  0.  (h6) 

By  comparison  with  (25)  we  see  that   (it-6)  implies  that 


and 


k  =  n^^(<^).  (U7) 


X  =  a  +h;[k^(a)][t-t^(0)|.  (1^) 


Here  we  have  used  the  fact  that  iii'(ik  )  =  h'(k  ).   ()f5)  and  (k8)   now  yield 

In  order  to  apply  (11.12)  we  require  the  second  derivatives, 

^10.   -=  -^"  f -*o] '   ^ka  =  -1  ^  ^'^^  '  ^aa   =  ^^^  •  (50) 

The  derivatives  are  evaluated  at  the  stationary  point,  0^,^  "  [^  (°)f°]-     We  set 

Then  (50)  and  (5I)  yield 
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det(^^j)  =  .[(l-h't;)2  +  (t-t^)h"  h  v] 
«=  -TCl-h'  f  )'^  +  (t-t  )h"  h  t'T 


(52) 


Again  we  have  used  the  fact  that  +h'(ik  )  =  h'(k  ). 

For  t  =  t  we  note  that  detfd)  .)  =  -1.  Therefore,  in  a  neighborhood  of 
t  =  t  (and  in  fact  until  the  ray  reaches  a  caustic),  sig((})  .)  =  0.   Therefore 
we  see  from  (ll . 12 )  that  u  ~  u .  +  u  ,  where 


ia(ih^,a) 


2ih7-  \^^-KK^'  -^  (*-*o)^;v;]"'^'  ^^  [^^[*^^;-^^o)  (*-*o)  ] j  > 

(53) 
x  =  a+h;(t-t^). 

We  see  that    (53)  and  (UO)  agree  eiiractly  if  we  set 

±a(±  h  ,0) 

%±  =  0'      ^±(^)   ■         2ih"  •  (54) 

o 

Here  h^  =  h^[k^(a)J,  and  a(h,a)   is  given  by  (hk) .      For  the  case  of  variable 

coefficients   (5^;)  supplies  the  values  of  s  ^  and  z^(cr)  necessary  to  conrplete  the 

construction  of  the  asyn^otic  solution.     In  that  case  we  must  set  h     »  h   fk  (o),a] 

o    o"-  o    •> 

in  (54). 

An  interesting  special  case  of  the  results  of  this  subsection  is  a 
source  function  of  the  form 

f(t,x;X)  ,  6[X(t.t^)]7(::)  =  ^  6[t.t J::)]7(x)  .  (55) 

In  this  case  g(t',a)  =  5(t' )7(cr)  and  a(h,cj)  =7  (a),  it  follows  that 
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z+(o)  =  ±?<<J)/2ih^        and 


(56) 


x=   aH.h;(t-tj. 


If  f  is  of  the  form  (55)  where  t  is  a  constarrt,  then  (25)  implies  that  k  (cr)  =  0. 
Hence  h  =  [b|,  x  =  o,  and  (56)  becomes 

u^-  ^  sin  [Xb(t-t^)].  (57) 

A  source  function  of  the  form  (21)  is  said  to  have  rapid  time  variation. 
One  of  the  form 

f(t,xiM  =g[t,x[x.xjt)]j  (58) 

is  said  to  have  rapid  space  variation.      Since  the  treatment  of  (58)  is  similar 
to  that  of  (21)  we  shall  not  examine  it  further.     A  special  case  of  (58)  is  the 
moving  jxjint  source, 

f(t,x;X)   =7(t)6|x(x.x^)]    =  ^  7  (t)s[x.x^(t )]  .  (59) 
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k.S     Oscillatory  point  source. 

A  source  function  of  the  form 

\^f(t,x;X)  =  x26(Xx)p(t)e"^^^^*^  =  X6 (x)p(t)e"^^'l^*^  (60) 

will  be  said  to  represent  an  oscillatory  point  source  located  at  the  origin. 

In  constructing  the  asymptotic  solution  of  (1,  2,   60)  it  is  natural  to 
assxime  that  waves  emanate  from  the  t-axls:   x  =  0,  t  =  t  (t  <  t)  with  a 
phase  given  at  x  =  0  by 

s(t,0)  =  s^(t)  =  -q(T).  (61) 

Differentiation  of  (6l)  with  respect  to  t  leads  to 


w  =  h(k)  =  +  h  (k)  =  q'(T);    h  =  sjc\^   +  b' 


(62) 


o 

This  equation  has  real  solutions  k  if  and  only  if  [q' (t)]   >  b  (O). 
For  definiteness  we  shall  assume  that* 

q*(T)  >  13(0).  (63) 

If  (63)  is  satisfied  then  (62)  has  real  solutions  k  only  for 


h  =  +  h^,  i6k) 


and  then 


(65) 


As  in  the  preceeding  section,  the  explicit  formula  for  k  (t)  is  easily  obtained 
but  will  not  be  used.   The  corresponding  values  of  g  at  the  source  are  clearly 
given  by 

g  =  h;  (+k^)  =±h'(k^).  (66) 

The  value  of  the  function  z  is  determined  by  its  value 


*In  terms  of  the  physical  interpretation  of  (l)  given  in  the  introduction,  the 
condition  (63)  is  that  the  frequency  of  the  source  should  be  greater  than  the 
plasma  frequency  ^.      It  is  well-known  that  waves  of  frequency  less  than  (j)  are 
exponentially  damped  in  a  plasma. 
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z(t,0)  =  z_|_  (t) 
at  the  source.   If  z   (t)  were  known  the  solution  could  be  completed  by  the 
methods  of  section  (2).   As  usual,  we  shall  determined  the  value  of  z   (t) 
by  computing  the  construction  for  the  case  of  constant  coefficients  and  by 
comparison  with  the  expansion  of  the  exact  solution. 

If  the  coefficients  are  constant,  the  rays  are  given  by 


X  =  +  h^[k^(T)](t-T),       t^  <  T  < 
and  are  illustrated  in  the  following  figure. 


(67) 


Those  corresponding  to  the  upper  (lower)  sign  in  (67)  proceed  to  the  right 
(left)  with  increasing  t.   As  usual,  s  is  given  along  each  ray  by 

s  =  s^(t)  +  t(t-T)  =  -q(T)  +  (t-T)[k^(T)h;(k^)  -  q'(T)J.  (68) 


(Here  we  have  made  use  of  the  equations  +h'(+k  )  =  h' (k  )  and  h  (k  )  =  q'(T).) 

^  ^  —  O  —  O        00  0^0      tL  \   /   / 

The  function  z  is  given  along  each  ray  by 

Z(t)  =  Z^(T)[j(T)/j(t)]^/2  (g^) 

where  j(t)  =  j(t,T)  =  dx/dcr.   In  order  to  calculate  j,  using  (67)  we  first 
obtain  the  value  of  k'(T).   This  is  easily  obtained  by  implicit  differentiation 


of  (62)  with  h  =  +h 


The  result  is 


k;(T)  =  q"(T)/h'(k^). 


(70) 
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j(t,T)  =  +[h'(k  )J-^fh"(k  )q"(T)(t-T)-[h'(k  )]2] 


and 


j(t,t)  =+h;(k^). 


{G^)    now  yields 

z(t)  =  z_^(t) 


(t-T)q"(T)h;;(k^) 


iK^K^y 


-1/2 


(71) 


(72) 


(73) 


hence  the  asymptotic  solution  of  the  intial-value  problem  (1,  2,    6o)  is  given 
parametrically  in  terms  of  t  "by 


u  ~  z_^(t)  1 


(t-T)q"(T)h;;(k^) 


[h:(k)]' 


■1/2 


JiXf-(l(T)+(t-T)[k^(T)h;(k^)-q'(T)]] 


X  =  +  h;[k^(T)J(t-T);  t^  <  T  <  t.  (7^+) 

The  exact  solution  of  the  problem  is  given  by  (1.23).   We  set 


rp(T),  t  <  T  <  t"1 

p(t)  =  >  ^      Y 

t,0   ,  otherwise  J 

Then  if  ve  insert  (60)  into  (1.23)  we  obtain 

h=+h  •     -       ^  ' 


(75) 


(76) 


;k,T)  =  -q(T)  +  kx  -  (t-T)h(k), 


(77) 
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We  now  apply  the  method  of  stationary  phase  to  (76).   Stationary  points  are 
determined  by  the  conditions 

(f)j^  =  X  -  (t-T)h'(k)  =0;         <j)^  =  -q'(T)  +  h(k)  =  0.  (78) 

Again  assuming  that  (63)  is  satisfied,  we  see  that  (78)  is  satisfied  only  for 
h  =  +h  and  then  k  =  +  k  (t).  Thus  only  one  term  of  (76)  yields  a  stationary 
point,  and  at  this  point 

X  =  (t-T)h;[+  k^(T)]  =  +  (t-T)h;[k^(T)].  (79) 

It  is  now  an  easy  matter  to  compute  the  matrix 

(i  )  =  (80) 

'  V.k     kJ 

at  the  stationary  point.   Then 

det(<^^.)  =  (t-T)q"(T)h^[ikJ-[h;(±k^)]2  =    (t-T)q"(T)h;;(k^)-[h;(k^)j2.  ( Qi) 

In  a  neighborhood  of  t  =  t  the  determination  is  negative,  hence  sig((j)^  .  )  =  0. 
The  stationary  phase  formula,  (11.12)  now  yields 

^-^  It  iH^\)f   -    (t-T)q"(T)h;;[kJ}-l/2  £(t1_  expjlx(-q(x)+(t-T)[k^(T)h'(k^)-q'(T)]]|: 

x  =  +  h^[k^(T)](t-T);      ^i  <  T  <  t.  (89) 

Comparison  of  (7^+)  and  (89)  shows  that  they  agree  exactly  if  we  set 

(^)  =   ish) .  (90) 

-       2ih'rk(T)]ci'(T) 
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For  the  case  of  variable  coefficients,  (90)  supplies  the  formula  for  z   (t) 
necessary  to  complete  the  construction  of  the  asymptotic  solution.   Important 
special  cases  of  the  results  of  this  section  are  obtained  by  setting 

q(t)  =  vt  [v  >  b(0)]  (91) 

and 

q(t)  =  vt   [v>b(0)];      P(t)  H  P^.  (92) 

(92)  yields  a  time-harmonic  point  source  of  frequency  W,  and  (91)  may  be 
interpreted  as  an  amplitude-modulated  time-harmonic  source. 
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k.  h    Slcfwly  varying  source 

As  a  final  example  of  our  method  of  solving  the  initial  value  problem  (l,2) 
we  consider  the  slowly  varying  source  function 

x2f(t,x;  X)  =  g(t,x).  (93) 

We  assume  that  g(T,x)  =  0  for  t^  <  t  where  tp  is  a  constant,  and  we  examine  the 
solution  u(t,x)  for  tp  <  t.  By  virture  of  (2)  we  may  define  gCx^x)  5  0  for 
T  <  t,  .  Thus 

g(t,x)  =  0  for  t  <  t^  or  t^  <  t.  (9k) 

Then,  in  general,  g(t,x)  will  have  a  discontinuity  on  the  line  t  =  t^,  and 
perhaps  on  other  curves  of  the  form 

t  =  t^(o),  X  =  a.  (95) 

As  in  section  (U.2)  we  take  t'(a)  >  0,  for  simplicity. 

We  assiune  that  each  discontinuity  curve  (95)  produces  tvro  waves  and  that 
on  the  curve  s  is  constant  for  each  wave.  We  shall  calculate  the  waves  produced 
by  one  discontinuity  curve  (9  )•  The  calculation  is  identical  to  that  of  section 
k.2   and  we  obtain  the  two  waves  u^  given  by  (UO).  Again  we  must  impose  the  condition 
(27)  and  z^(ct),  s  ^  are  undetermined.  In  order  to  determine  them  we  use  the  exact 
solution  (1.23).  By  virtue  of  (qH)   we  may  write  it  in  the  form 

U  .  r-    I^I  /  dl./  d.  /  dx  SLL^      UQ,(x-C).(t.T)h]^  ^  ,  ^^(^,     f^g, 

' '•    ■■  —00       —00         .00 
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Now  let 

7{o)   =  [g]    =g[t^(o)-K),    a].g[l^(a).0,aj  (97) 

be  the  magnitude  of  the  jump  discontinuity  of  g.     Then  by  integration  by  parts 

f     ^L^  e^^^^  dT  ^-    r(a)e^'*o("^^   .  (98) 

-00 

Next  we  insert   (98)  in  (96),     The  res\ilt  is 

00  OJ 

u~^      -^/     ^X^^'     7(a)e^'^^^>^^  (99) 


where 


^(k,o)  =k(x-a)  .  h|t.t^(a)].  (100) 

The  integral  (99)  nay  be  evaluated  asymptotically  by  the  method  of 

stationary  phase.  Since  (90)  is  of  the  same  from  as  (^3)>  we  may  obtain  the 

afh  a)    ^{^   >a) 

result  at  once  by  replacing  ^^  j*  '  =    rt-2 by  -7(«^)  in  (53).  We  may  then 

o 
determine  c  .  and  z^(cj)  by  comparison  with  (hO).     In  this  way  we  obtain 

s^±  =  0,  z^(a)  =.jl£2  =  .|j;g],  (101) 


and  these  results  may  be  used  to  obtain  the  asyinptotic  solution  of  the  initial 
value  problem  with  variable  coefficients. 
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Appendix  I.   Exact  solution  of  the  initial  value  problem  for  constant  coefficients 
We  first  consider  the  initial  value  problem  for  the  homogeneous  equation, 

c\  -  u^.  -  xVu  =  0,  (1) 

XX    tt  ' 

u(0,x)  =  u^(x;  X),   u^(0,x)  =  u^(x;X).  (2) 

We  consider  a  solution  of  (l)  in  the  form  of  a  superposition  of  plane  waves, 

(3) 

r  iX(kx±h  t)  /op? 

u(t,x)   =  y~     /       a^(k,X)e  "^    dk;     h^  >   >/=  k  "^^   • 

In  order  to  satisfy  the  initial  conditions  (2)  we  set 

u(0,x)  =  y /a^''^^  dk  =  u^(x;X), 

u^(0,x)   =  Y'J    >iXh^a/^^^  dk  =  U3^(x;X). 
Then,  by  means  of  the  fourier  transform  formulas, 
g(x,X)  =    jra(k,X)e^^^       dk;     a(k,X)   =^    /g(x,X)e-^^^  dx, 
we  may  inVert   (h)  and  (5).     This  yields 

a^+a_   =^    /u^(x;  X)e-i'^^  dx,  (?) 

k  /u,(x;x)e-^^^dx,  (8) 


(5) 


•+  ■".         2i:h 
o 


and  therefore 


(6) 


a,(k,X)  =  I-     /KuJx;X)±  t-  u^(x;X)]e-^^^^  dx.  (9) 

o 
Finally  the  solution  of  (1,2)  is  obtained  by  Inserting  (9)  in  (3). 
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This  yields 


u(t,x)  =l-y      I     6k   I     dafxu  (a;X)±  J-   u(a;X)]e  L       oj 
L-      -co       .0=  I  o(k)       J 

(10) 

The  solvttion  of  the  initial  value  problem  for  the  inhomogeneous  equation 

can  be  reduced  to  the  solution  of  (1,2)  and  the  solution  of  the  problem  for  the 

inhomogeneous  equation  vith  zero  initial  data.  Therefore  ^  e  consider  the  problem. 


u(t^,x)  =  u^(tpx)  =  0.  (12) 


Here  t.  is  a  constant.  In  order  to  solve  (11,12)  ve  set 

u(t,x)  =  fa.(t,k)e^^^''  6k}   f(t,xjX)  =  f^(t ,k}\)e^'^^''   dk.      (I3) 

Then  (ll)  yields 

a^^  +  >^V  a  =  -X^<|,,  Cl'^) 

and  from  (12), 

a(t^,k)  =  a^(t^,k)  =  0.  (15) 

To  solve  the  initial  value  problem  (lU,15)  for  the  inhomogeneous  linear  ordinary 
differential  equation  (lU)  ■  e  apply  the  method  of  "variation  of  parameters"  to 
the  tvo  linearly  independent  solutions  exp  [  ±i^h  t  |  of  the  homogeneous  equation; 
i.e.  ve  set 


iXh  t  ^   -^^V  ,,.. 

a  =  v^e   o  +  Vpe  (16) 


and  insert  (l6)  in  (ik).      In  so  doing  "e  require  that 


iXh  t      -iXh  t 


and  ve  obtain 

iXh  t         -iXh  t     p 
i^Vlt^   °  -i^V2t^    °  "-^^  (1^) 

(17)  and  (18)  may  no'  be  solved  for  v^  and  v^.  and  then  integrated  to  obtain 

,2     p   -iXh  s  \-        r      i^h  s 

°  t3^  °h 

We  no^-  insert  (19)  in  (16).   This  yields 

.     P        liXh  (t-s) 
a»y~  ±  2j^   /   e    °      .f)(s,k;^)ds,  (20) 

and  it  is  easily  seen  that  (15)  is  satisfied.  By  inverting  the  fourier  transform 
(13)  ire   see  that 

H^,^;^)   =  |j   rf(s,y;X)e-^^^y  dy.  (21) 

,5,2    P     p  iXrky±(t-s)hri 

^^T-'&  y  ^^1      ^^^  ^-'fCs^yjX)  (22) 


Then 
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and  by  insertinc  ('52)  in  (I3)  "e  obtain  the  solution  of  (11,12), 

,(,„.x,  =y-  ^    f   ^f   d,.    r\3.C(^.>:i^->-e'^&(-^'-(*-='^]-     (23) 

^ —  loo       -00        t^ 

Here  h(k)  =  fl^,(k). 
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Appendix  II.   The  method  of  stationary  phase  in  n  dimensions 

In  this  appendix  "e  derive  a  simple  fonmila  for  the  leading  term  of  the 
asymptotic  expansion  (for  X  -»»)  of  a  function  f(X)  defined  by  a  certain  type 
of  n-dimensional  integral.   The  method  'e  use  is  the  method  of  stationary  phase. 
Although  numerous  references  for  this  method  exist,  the  n-dimensional  formula 
ve  require  is  not  readily  available.  Ve   -ill  derive  it  by  a  method  ■  hich  is, 
if  not  entirely  rigorous,  relatively  simple.  Let 

f(X).  r.(t, Ve"*''^'-'*"'dt,...dt„  =  f.(T)e^H(TV.        („ 

A  stationary  point  is  a  point  T  =  S  =  (s^, . . ,  s  )  for  -.rhich  the  real- valued  phase 

function  ^(T  )  Is  stationary,  i.e.  i  (s)  =  ^  (s)  =  Oj  v  =  l,...,n.  We  assume 

V       ^ 
that  the  domain  of  integration  D  in  euclidean  n- space  E  contains  exactly  one 

stationary  point,  and  that  this  point  lies  in  the  interior  of  D.   (if  D  contains 

several  stationary  points  it  can  be  subdivided  into  sub-domains  each  containing  one.). 

A  Taylor  expansion  of  ^(^H)   about  the  stationary  point  then  yields 

KT)  =KS)  +|(t^-s^)(t^-s.).f)^^  +o[|T  -SpJ.  (2) 

Here  ^  ^  =  <|>j.  ^  (S),  and  the  summation  convention  •"ith  respect  to  repeated  indices 

"  J 
from  1  to  n  is  employed. 


We  now  introduce  a  translation  and  rotation  of  T-space  so  that  the  point 
T  =  S  is  the  origin  and  the  matrix  i^y  .)    is  diagonal.   Then 

6    .   =  T     &    ..  (3) 

Tvj    V  vj  ^^J 
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Here  r  ,...r  are  the  (real)  eigenvalues  of  ((j)^.).   It  is  assumed  that  (<|)  .)  is 

1/2 
non-singular.   Hence  no  eigenvalues  vanish.   We  nov  set  y  =  \  '      and  introduce 

the  transformation 

Then 

.(t^-s,)(t.-s.)^^.  =  7\t.r,^,.   =   x^x.r^S^.  =   x^r^  (5) 

and  from  (2) 

exp[iX^(T  )]  =  exp  {iX<f,(s)  +  ^^y^   +  0(X-^/2)  ]   =  exp  [iX^(s)  +  ^r^^x^  ^  [l-K)(X-V2)l 

(6) 
Furthermore  g(  T  )  has  the  expansion 

g(  T  )  =  g(s)  +  0(X-V2)  ^  g(3)ji+  o(X-V2)J  ^^j 


i  =  7"  =  X-/2  (8) 

^•Te  no-^  insert  (6),    (7)  and  (8)  in  (l).   This  yields 

f(X)  ~  gC3)e^H(S)f.  .   ^.  ^^-n/2j'h  \     ^^  ^^^ 

Here  D'  is  the  image  of  D  under  (k).     Ve  assume  that  a  negligible  error  is 
introduced  if  ve  replace  D'  by  the  infinite  X- space.  Then 


k=l  ^00 
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Ve   now  apply  to  (lO)  the  standard  integral  formula 


Here  sgn  r  =  ±  1  if  r  ^  0.   Then  (9)»ClO)>  and  (ll)  yield  the  required  formula 
f(M~(^  [l  det(<>^j)|]-^/^  g(S)e  "^   }   X  -.».        (12) 

In  (12),  sig  C<t'y-)  denotes  the  "signatxire"  of  the  matrix  (^    .)  defined  by 

n 
sig  ((S^j)  «  y^  sgn  r^,  (13) 

k=l 
vhere  r^,...,  r  are  the  eigenvalues  of  ((}>  .).  We  have  also  used  the  identity, 

We  note  that  both  |det(<))^.)|  and  sig(<f)^.)  are  invarient  under  translation 
and  rotation  of  the  co-ordinate  system,  hence  the  result  (12)  is  valid  for  the 
original  co-ordinates. 
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